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VANISHING VISCOSITY AND SURFACE TENSION LIMITS OF 
INCOMPRESSIBLE VISCOUS SURFACE WAVES 

YANJIN WANG AND ZHOUPING XIN 


Abstract. We consider the dynamics of a layer of viscous incompressible fluid, subjected to 
gravity and surface tension forces. The upper boundary is a free boundary, and the lower 
boundary is a fixed boundary on which the Navier-slip boundary condition is imposed. It is 
proved that there is a uniform time interval on which the uniform estimates independent of 
both viscosity and surface tension coefficients can be established. These allow one to justify 
the vanishing viscosity and surface tension limits by the strong compactness argument. In 
particular, a unified local well-posedness of the free-surface incompressible Euler equations 
with or without surface tension can be obtained by the inviscid limits. 


1. Introduction 


1.1. Formulation. We consider the motion of an incompressible viscous fluid, subject to the 
influence of gravity and surface tension forces, in a moving domain 

Q(t) = {x € R 3 | — b < x 3 < h(t,x 1,2:2)} • (1.1) 

The lower boundary of Q(t) is assumed to be rigid and given with the constant b > 0, but the 
upper boundary is a free surface that is the graph of the unknown function h : R + xR 2 —>• R. The 
fluid is described by its velocity and pressure, which are given for each t > 0 by u(t, •) : Q(t) —>• R 3 
and p(t, •) : 12(f) —>• R, respectively. For each t > 0, ( u,p , h) solve the following problem for the 
incompressible Navier-Stokes equations 


dtu + u ■ Vu + Vp — eAu = 0 
V • u = 0 

< pn — ‘2sSun = ghn — aHn 
dth = u ■ N 

u 3 = 0, (Su(-e 3 ))i = -KUi, i = 1,2 


in 12(f) 
in 12(f) 

on {z 3 = h(t,xi,x 2 )} 
on {x 3 = h(t,xi,x 2 )} 
on {x 3 = -b} 


( 1 . 2 ) 


for Su = -7 (Vu + Vu 4 ) the symmetric part of the gradient of u and n = N/|N| the outward 
unit normal of the free surface with N = {—d\h,—d 2 h,l) t . e > 0 is the viscosity, g > 0 is 
the strength of gravity, cr > 0 is the surface tension coefficient and k is the friction coefficient. 
Finally, H is twice the mean curvature of the free surface given by the formula 


H = V- 


Vh 


A + lv/il 



(1.3) 


The kinematic boundary condition, the fourth equation in 111.21) . implies that the free surface 
is adverted with the fluid, and the dynamic boundary condition, the third equation, states 
the balance of stress tensor on the free surface. We have imposed the Navier slip boundary 
condition on the lower fixed boundary. Note that in (11.21) we have shifted the gravitational 
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forcing to the boundary and eliminated the constant atmospheric pressure, p a tm, in the usual 
way by adjusting the actual pressure p according to p = p + gx 3 — p a tm- 

The initial surface is given by the graph of the function h{ 0) = ho : M 2 —>• R, which yields 
the initial domain H(0) on which the initial velocity u(0) = uq : 0(0) —>• M 3 is specified. It will 
be assumed that ho > —b and that (uo,ho) satisfy certain compatibility conditions. 

The movement of the free boundary and the subsequent change of the domain create numerous 
mathematical difficulties. To circumvent these, as usual, we will transform the free boundary 
problem under consideration to a problem with a fixed domain and fixed boundary. Consider a 
family of diffeomorphism <f>(t, •) of the form 


$(t,-) : 0 = R 2 x (-6,0)^ 0(f) 

(y,z) ->• {y,v(t,y,z))- 


(1.4) 


tp is chosen so that d z <p > 0 which ensures that <h(f, •) is a diffeomorphism. For the fluid domain 
under consideration, p can be chosen as 

<fi(t,y,z) = z + rj{t,y,z), (1.5) 

where y is a chosen extension of h onto {z < 0 } defined by 

V(t, z) = (l + exp (A\£\z)h(t, £). ( 1 . 6 ) 

Here : stands for the horizontal Fourier transform with respect to the y variable. It is verified 
in Proposition IB. II that for given ho if the number A > 0 is chosen sufficiently small, then 


d z p(0, y, z) > Co > 0 in O. 

Then one can reduce the problem into the fixed domain 0 by setting 

v(t, y, z ) = u(t, y, z)), q(t, y, z) = p(t, ${t, y, z)) in H. 


Set 


df = di- jr^-d z , i = t, 1 , 2 , d% = = ~^—d z 

d z p 6 d z jp 


such that 

diU o $(t, •) = dfv, i = t, 1,2,3. 
Then by the change of coordinates USD, the problem (11.21) becomes 


dfv + v 

■ + V^q 

- eA^v = 0 

in 

n 

V 1 * 3 -v = 

0 


in 

n 

qn — 2eS‘ p vn = ghn - 

- aHn 

on 

{z 

d t h = v 

• N 


on 

{z 

V3 = 0 , 

(S^ve^i = 

KVi, i = 1,2 

on 

{z 

(v, h ) \t= 

0 = (vo, ho). 





Here we have naturally written (V^ = df, = df df , • v = dfv, and S v v 

Note that • S^v = for vector fields satisfying • v = 0. 


(1.7) 

( 1 . 8 ) 


(1.9) 


5(V^ + 


1.2. Previous works. Free boundary problems in fluid mechanics have been studied inten¬ 
sively in the mathematical community. There are a huge amount of mathematical works, and 
we only mention briefly some of them below. We may refer to the references cited in these 
works for more proper survey of the literature. For the incompressible Navier-Stokes equations, 
we refer to, for instance, Beale [ 6 ], Hataya [22]) Guo and Tice HSJIIBEU for the well-posedness 
without surface tension, and Beale [7], Tani [39], Tanaka and Tani [50] for the well-posedness 
with surface tension. Those well-posedness results are strongly based on the regularizing effect of 
the viscosity, and the solutions are shown to be global for the small initial data [ 22 ] 1201 EH 0 gO]. 
Note that the surface tension only has a regularizing effect on the free surface, and its effect 
serves to enhance the decay rate, see [20 1 [21 ] for more discussions. For the incompressible 
Euler equations, the problem becomes much more difficult. The early works were focused on 
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the irrotational fluids, which began with the pioneering work of Nalimov [31] of the local well- 
posedness without surface tension for the small initial data and was generalized to the general 
initial data by the breakthrough of Wu 0H H2 for the case without surface tension and by 
Beyer and Gunther [9] for the case with surface tension. For the irrotational inviscid fluids, 
certain dispersive effects can be used to establish the global well-posedness for the small initial 
data; we refer to Wu 03JIS], Germain, Masmoudi and Shatah m, Ionescu and Pusateri [25] 
and Alazard and Delort [Ij for the case without surface tension, and Germain, Masmoudi and 
Shatah m and Ionescu and Pusateri [26] for the case with surface tension. For the general 
incompressible Euler equations without the irrotational assumption, only local well-posedness 
results could be found. The first local well-posedness in 3D was obtained by Lindblad [29] 
for the case without surface tension and by Coutand and Shkoller HD for the case with (and 
without) surface tension. We also refer to the results of Shatah and Zeng [36] and Zhang and 
Zhang 07] . 

Various approaches are used to prove those well-posedness results mentioned above, depend¬ 
ing on whether viscosity or surface tension is presented or not. It is then very natural and 
interesting to study the asymptotic behavior of vanishing these two parameters in the equa¬ 
tions. The vanishing viscosity limit for the Navier-Stokes equations is a classical issue. When 
there is no boundary, the problem has been well studied; we refer to Swann EZ], Kato EZ], 
DiPerna and Majda Ezma, Constantin [10] and Masmoudi [51j for example. However, in the 
presence of boundaries, the situation is more complicated and the problem becomes challenging 
due to the possible formation of boundary layers. In a fixed domain with the no-slip boundary 
condition, there is formation of boundary layers in the vicinity of the boundary and the solu¬ 
tion u £ of the Navier-Stokes equations is expected to behavior like u £ ~ u° + U (t, y, z / yfe) (we 
assume the boundary is locally given by z = 0), where vP is the solution of the Euler equations 
satisfying only the impermeable boundary condition and U is some profile. In view of this small 
scale behavior, it is impossible in general to get uniform strong estimates in any Sobolev spaces 
containing normal derivatives. Consequently, the vanishing viscosity problem with the no-slip 
boundary condition is widely open except the works by Asano [5] and Sammartino and Caflisch 
m in the framework of analytic initial data and and the work by Maekawa m for the initial 
vorticity located away from the boundary and a recent work by Guo and Nguyen [18] for a 
steady flow over a moving plane. However, when the no-slip boundary condition is replaced by 
the Navier slip boundary condition, the situation becomes better. Indeed, now the solution is 
expected to behavior like u £ ~ u° + y/sU(t, y, z/y/e)] the amplitude of the boundary layer is 
weaker, and one can hope to get an uniform estimates involving one normal derivative. In this 
case, the vanishing viscosity limit in 3D has been justified rigorously in Iftimie and Planas [23], 
Iftimie and Sueur [23], Masmoudi and Rousset [32] and Xiao and Xin 06]. Furthermore, for 
some special types of Navier boundary conditions or boundaires, uniform estimates in higher 
order Sobolev spaces can be obtained, see Xiao and Xin 05] and Beirao da Veiga and Crispo 

El- 

Going back to the free-surface incompressible Navier-Stokes equations, since the dynamic 
boundary condition can be viewed as the same type of slip boundary conditions, one has the 
hope to establish the vanishing viscosity limit. For the case without surface tension and there 
is no boundary below the fluid, Masmoudi and Rousset [33] justified the inviscid limit by using 
the framework of their earlier work [32] and some additional techniques. Later, Elgindi and Lee 
[13 ] discussed the same problem for the case with surface tension, however, some key points in 
their arguments are not clear to us. On the other hand, for the free boundary problems, it is 
also interesting to show the vanishing surface tension limit. This is supposed to be somewhat 
simpler than the inviscid limit problem since the equation on the free boundary is defined 
without boundary. Yet, one needs to develop the well-posedness which is uniform with respect 
to surface tension, and generally this is nontrivial. These have been done for the irrotational 
Euler equations, see and references therein; for the general Euler equations, a priori 

uniform estimates have been derived in [39] i while for the Navier-Stokes equations with small 
initial data, we refer to the recent work of Tan and Wang [38] . The purpose of this paper is 
to derive the uniform estimates of solutions of the system (11.21) (equivalently, (11.91) 1 on a time 
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interval independent of both viscosity and surface tension coefficients. These allow us to justify 
the vanishing viscosity and surface tension limits by the strong compactness argument. As a 
byproduct, we can get a unified local well-posedness of the incompressible free-surface Euler 
equations with or without surface tension by the inviscid limits. 


2. Main results 


2.1. Statement of the results. We shall use Sobolev conormal spaces on as [32l 1331 . Set 


Zi = di, i = 1, 2, Z 3 = z(z + b)i%, 

which are tangent to <9fL The Sobolev conormal space H k a is defined as 

Him = {/ € L 2 m, z a f g L 2 m ,« g n 3 , m < k }, 

where Z a = Z^ 1 Z^ 2 Z^ 3 , with norm dehned as 

ii/iu = E ii z “/ib. n/ii = ii/ii 0 = ii/ib. 

|ck| <fc 

Similarly, Wco°° is defined as 

W*-°°(fi) = {/ G L°°m, Z a f G L°°(Q), a G N 3 , \a\ < k} 

with norm 

ii/iu,«, = E ii z ”/iu~ • 

|O'| <k 

H k and W k, °° will denote for the usual Sobolev spaces on ST, and | • | s and | • | S!00 stand for the 
standard Sobolev norms on R 2 . To estimate the time derivatives of the solution, as m , we 
introduce also the spatial-time Sobolev conormal norms on S2 as: 


m 

= E I®!/ 

£=0 


2 

m+k—i 


i,o , and 


m 

E 

e=o 


Wtf 


m—i. oo 


( 2 . 1 ) 


and the spatial-time Sobolev norms on R 2 : 


,’ i/ix™ = i/ix m ’°. and i/Iy™ = 1 5 */ 

L —' 1 m+s—t *— 


£=0 


£=0 


m—t,oo 


( 2 . 2 ) 


Note that in these definitions k and m are assumed to be non-negative integers, but s is allowed 
to be any real number, typically, halfs. In the following we will abuse the notation ^ for [?(■] • 
The aim of this paper is to get a local well-posedness result for strong solutions of (USD in 
Sobolev conormal spaces which is valid on an interval of time independent of e, cr G (0,1]. Note 
that such a result will also imply the local well-posedness for the Euler equation with or without 
surface tension. As it is well-known that when there is no surface tension a Taylor sign condition 
on the free boundary is needed to get local well-posedness for the Euler equation; when there 
is surface tension, no such condition is needed. By the change of coordinates (11.41) . the Taylor 
sign condition reads as 

— d%q + g > Co > 0 on {z = 0}. (2.3) 

In the below, N and n are extended to fl by 

N(«,y,z) = (-Oi</?(t,?/,z),-« 92 </?(t, 2 /, 2 ), 1 Y and n = N/|N|. (2.4) 


Note that N(f,y,0) is indeed the outward normal to the free surface defined before, and 
N (t,y,—b) = e 3 . Define II = I — n ® n, and let x de a smooth function which takes the 
value zero in the vicinity of {z = 0} and one in the vicinity of {z = —b}. Let ||-||^px be the 
norm of the space L p ([ 0, T]; X). Then the main result in this paper is stated as follows. 
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Theorem 2.1. Let m > 14. Assume that the initial data (/iq a , v E q U ) is given such that 

\h £ ’*mlm +a +e\h^m 2 x ^ (2.5) 

+ 11^(0)11^ + ||a 2 p^(0)||2 m _ 1 + 1 | 4 ^( 0)||^ +2 +e\\d zz v^ml~ < Ro, 

and that initially the compatibility condition H £,rJ v £,a n £,cr — Kxv e,a ) = 0 on {z = 0,-6}, 

the diffeomorphism m and the Taylor sign condition (IMD on {z = 0} hold. There exists 
T > 0 and C > 0 such that for every e,cr € (0,1], the unique solution (v £ ’ a , h £ ’ a ) of (11.91) 
satisfies the estimate: 


—M + a L5?X™-L 2 + £ 


( 2 . 6 ) 


+ ll t ’ £ ’ IIl^X ™- 1 ' 1 + II^W^’ llz,^?X m_2 "b IIj^y^T " 1 " 2 £ ^® zzvE ’ IIL^?L°° 

+ 0-2 l^ £ ’ ff | 2 ^ x m-l,| + 6 IIV^Ili^-M + £ l|V<9 2 P e,CT ||^ X m-2 < C. 


It should be remarked that the reason for the L-in-time estimates of ||5 2 u £: ’ f 7 (t)|| Xm _i rather 
than in L°° stated in Theorem 12.11 is related to the boundary control of the vorticity for 
the Navier-Stokes system (be., viscous boundary layers). Also, owing to the less regular¬ 
ity of the pressure as will be shown in our proof later, which is due to the presence of sur¬ 
face tension, we can only prove the L-in-time estimates of the highest time derivative term 
\df n h £ ’ a (t )\ L 2 + y/o \df 1 h £,ar (t)\ H i + \fe\df l h £,C7 (t)\ i + \\df 1 v £,a (t)\\ L 2 . As immediate conse- 
quences of the uniform estimates of Theorem 12.11 one establishes easily, by standard compact¬ 
ness arguments, the justification of the inviscid limit, vanishing surface tension limit, and any 
their combinations. In particular, we can obtain a local well-posedness of the free surface Euler 
equations, independent of surface tension, by the inviscid limit. Furthermore, for the Euler 
equations, by using the equation for the vorticity one can improve those L 4 -in-time estimates 
to be in L°° and also recover the standard Sobolev regularity. 


2.2. Strategy of the proof. The main step in the proof of Theorem 12.11 is to derive the a 
priori uniform estimates on a time interval small but independent of e, a € ( 0 , 1 ] for a sufficiently 
smooth solution of the equations (11.91) . Our approach is strongly motivated by the strategy of 
Masmoudi and Rousset [33] where the vanishing viscosity limit was justified for the problem 
without surface tension. However, there are several new difficulties arising in the presence of 
surface tension. Indeed, it has been already known for the free-surface incompressible Euler 
equations that the problem with surface tension is more difficult than the problem without 
surface tension in certain sense, see Coutand and Shkoller El for some discussions. This is 
even so when one wants to show the inviscid limit for the free-surface incompressible Navier- 
Stokes equations. When showing the inviscid limit for each fixed a > 0, one finds that surface 
tension only serves to provide the improved regularity of the free surface. Indeed, the presence 
of surface tension makes the problem more difficult: first, the required nonlinear estimates are 
more difficult to close due to the mean curvature term; second, the less regularity of the pressure 
makes the arguments much more involved. Note that the less regularity of the free surface for 
the problem without surface tension was get around in [33] by using Alinhac good unknowns 
[2]. In this paper we will still use Alinhac good unknowns also without taking advantage of the 
improved regularity of the free surface provided by any fixed surface tension; nevertheless, we 
will develop some ideas to overcome the difficulties caused by the presence of surface tension as 
illustrated below. As a consequence, we will be able to deal with both the vanishing viscosity 
and surface tension limits. 

For notational convenience we shall suppress the subscripts e and a below. Let J\f(T) be 
the quantity appearing in the left hand side of ([2.61) . while Q(T ) be the first two lines in ([2.61) . 
The crucial point is thus to get the closed a priori estimates of J\T(T) on a small time interval 
independent of £ and a in terms of the initial data. We start with the basic physical energy 
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identity: 

“ H 2 dV t + J z _ o 9\h\ 2 + 2a(^Jl + \V y h\ 2 - l) dy') (2.7) 

+ 2 s f {S^vfdVt + ^Ke f \v\ 2 dy = 0. 

J J z=—b 

Here dVt stands for the volume element induced by the change of variable (j 1.5 D : dVt = d z <pdydz. 

To get the estimates for higher order conormal estimates, one then applies Z a for 1 <\a\ <m 
to the equations (11.91) . Since the operators df involve V</?, the estimate of the commutator 
between Z a and df needs a control of ||Z Q V</?|| < |Z a /i| 1 / 2 . In the absence of the surface 
tension this yields a loss of 1/2 derivative, and the difficulty was overcome in Masmoudi and 
Rousset m by using a crucial cancellation observed by Alinhac [2]. The idea is to use the 
good unknowns V a = Z a v — dfvZ a rj and Q a = Z a q — dfqZ a r], and some cancellation occurs 
when considering the equations for V a and Q a which allows one to derive an L 2 type energy 
estimates similar to (12.71) : 

~ (7 m 2 dVt + f _ (g - dfq)\z a h\ 2 + a \S7 y Z a h\ 2 dy) +2 e f \S*V a \ 2 dV t (2.8) 

\t/ J 2—0 J J 

= - f aZ a H Cffz a 'n-Z a ~ a 'vdy- f C a (d)Z a qdV t + Y +J2 > 

J z= 0 i /i i J C2 0 cr 

1 “ 1 = 1 

where, using the symmetric commutator notation defined by (IA.3I) . 

d z ipC a {d) = [Z a , N, -d z v\ + [Z a ,d z rj, d x v x + d 2 v 2 ] . (2.9) 

Here denotes the terms that can be controlled in a similar way as the case without surface 

tension [33], and stands for the terms related to surface tension, after some delicate argu¬ 

ments, that can be controlled well with the improved energy estimates a \Z a h\f. The first two 
terms in the right hand side of (12.81) are singled out in order to indicate the main difficulties 
for the case with surface tension. Note that the regularity a \ Z a h\l in the energy is not enough 
to control the first term. Indeed, one needs a \Z a h\ 3 / 2 , *• e -t there is a loss of 1/2 derivative. To 
improve the regularity of h, one then resorts to the “dissipation” estimate by using the dynamic 
boundary condition: —oH = q — gh — 2eS lf vn • n on {z = 0}, which requires a control of 
\Z a q\_ 1 j 2 - On the other hand, it seems that the only way to control the pressure q is through 
the elliptic equation with the Neumann boundary conditions 

' A^q = -VV ■ (v • Wv) in Q 

< V^q ■ N = —dtv • N — (v y • V y )v • N + eA^v • N on {z = 0} (2.10) 

VA? • N = eA*v • N on {z = -b}. 

It should be noted that on the boundary {z = 0} one could not use the Dirichlet boundary 
condition q = —aH + gh + 2eS‘ p vn • n since there is no control for —aH yet, which is in 
contrast to the case a = 0 where the Dirichlet boundary condition is successfully used in [33]. 
Note that the appearance of dtv in (12.101) forces one to include the time derivatives in Z a . 
Then elliptic estimates for (12.101) provide a control of ||V(?|| xm -i, which then yields a bound for 
a l^l x m-i,§ • This implies in particular that we need to separate the estimates of (12.81) into two 
cases: «o < m — 1 and ao = m. For the former case, the main conclusion is that 

II^WIIX™- 1 . 1 + + CT \h(t)\xm-l,2 + £ |^(i)| 2 m _l,3 ( 2 -ll) 

< Co + A (Q(T)) (t + J* \d?h\\ + a\d?h\\ + \\d?vf + • 

The crucial step is to treat ([2.8D for the case Qo = m, which can be rewritten as 

( [ l ym | 2dV t+ [_ (g - dtq)\d?h\ 2 + a\Vyd?h\ 2 dy \ +2 e f \S*V m \ 2 dV t 

J Z —0 J J 12 


( 2 . 12 ) 
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/ ■ 


adVHmd t N ■ fj^-'vdy - / md t N • 


m— 1 


d'var<,dyiz+•£+•£ +y 


Here stands for the terms related to q, which are easier to estimate after some delicate 
arguments than the second term on the right hand side above. The main difficulties now are 
that there are no any estimates for d™q so that it is difficult to bound the second term, and hence 
that the regularity in the energy can no longer be improved through the dissipation 

so that it is difficult to bound the first term. Our approach to overcome these difficulties is to 
integrate by parts over t and z in an appropriate order to obtain a crucial cancelation. More 
precisely, we integrate by parts in z first and then in t the second term to obtain 


— f m<9fN • d™ l d z vd] n qdydz 

Jn 


(2.13) 


[ md t lSi ■d™- 1 vd™qdy+ [ md z (d,Nd™q) ■ d^vdydz 
J z =0 J 

f mdt N • d^-'vd^q dy + ^ f md z (d t Nd™- l q) ■ d™~'vdydz + V . 
Jz =o dt Jn ^o 


Note that one can not integrate by parts in t for the first term in the last line of (12. 13|) since we 
can not control d^v on {z = 0}. The crucial observation is that there is a cancelation between 
this term and the first term in the right hand side of (12.121) . Indeed, by the kinematic boundary 
condition, it holds that 


/ ■ 


ad^HmdtN -d^vdy- 


md t N ■ d^vd^qdy 


(2.14) 


'2 = 0 


' 2=0 


mdfN • d r t n -'v {gd^h + 2 sdj™ (S^v n • n )) dy = ■ 


The last difficulty in closing the conormal estimate for the case cho = m is due to the second 
term in the last line of (12.131) : the control of ||V(/|| Xm _i (especially involves the dissipation) is 
not pointwise in time but only L 2 . To get the estimates, we are forced to integrate in time twice 
and the main conclusion is that 


(lIM 2 + l^r^lo + * ml + e |^||) 2 ) 2 (2.15) 

< Co + A(Q(T)) + £ \d?h\l + a 1^1? + \\drv\\ 2 + \\d x v\\b»-i + * l^l|) ' ■ 


The next step is to derive the conormal estimates for d z v. As in [33], the main idea is to first 
introduce the equivalent quantity S n = n (S^vn — kxv), which satisfies a convention-diffusion 
type equation with the homogeneous Dirichlet boundary condition. The main difficulty is that 
Q(T) controls only y/e ||3 zz v|| i0 o. Since we have included the time derivatives, so if one followed 
the arguments of [ 33 ], which only involves the spatial derivatives, to estimate the commutator 
resulting from the viscosity term, one would need to control y/£\\d zz v\\ Y k for some high k > 1 . 
Our key observation here is that since in the vicinity of the boundary the solution behaves like 
v(t,x) ~ v°(t,x) + y/eU(t,y,z/y/F), so it is possible to control ed zz v (and even ed zzz v !). This 
can be done easily by using the first equation in (11.91) . We can then perform the L 2 type energy 
estimates to conclude that 

\\d z v(t)\\l m - 2 <C 0 + A(Q(T)) (t + j • (2.16) 

Note that the m — 2 order estimate above cannot be improved to be m — 1 due to the appearance 
of (V^) 2 ^ in the source term in the equation for Sn, which is only in X m ~ 2 . To get a better 
estimate, following [3311321115], one would proceed with the vorticity u = x v instead of 
S n . u) again satisfies a convention-diffusion type equation, but the main difficulty is that it 
does not vanish on the boundary and the boundary value is at a low regularity. To split the 
difficulty, for |a| = m — 1 we set Z a oj = uj^ h + oj where u:® h satisfies the nonhomogeneous 
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equation with the homogeneous boundary condition which can be handled by performing the 
L 2 type energy estimates and co^ h satisfies the homogeneous equation with the nonhomogeneous 
boundary condition. Note that u)® h solves exactly the same problem as [33] so that one can apply 
their estimates directly, though, the proof is quite involved. Combining these two estimates, 
one concludes that 

(^ll^)llx-i) 2 <C 0 + A(Q(T)) (t + J* lIMjU-.) • (2-17) 

The last step is to derive the L°° estimates of d z v and y/ed zz v. In this step, one needs to 
estimate a low number of derivatives of v, say y + 2, while the boundary is H m with m being 
as large as needed, it is convenient to use a normal geodesic coordinate system in the vicinity 
of the boundary so that the Laplacian has the simplest expression. Note that \fed zz v can be 
controlled in the same way as [33] . however, for d z v we will employ a different argument. Again, 
it is more convenient to estimate the equivalent quantity S n . After some computations, we find 
that p, an equivalent quantity of S n in the new coordinates, solves 

d t p + w -V p - ed zz p = 'H (2.18) 

for some source term TL and a vector field w with W 3 = 0 on the boundary. The main difficulty in 
the analysis of [S3] is a commutator between Z a and ed zz which is difficult to control when ap¬ 
plying the maximum principle. The main idea in [33] is to rewrite (12.181) into a one-dimensional 
Fokker-Planck type equation and then use the explicit representation of the solution. Our ob¬ 
servation here is that since we have included the time derivatives in Z“, it is easy to estimate 
the commutator by, roughly speaking, replacing the ed zz with dt by using the equations. Main 
conclusion in this step is that 

\\d z v(t)\\ 2 Y rn +2 + e ||3„u(i)||£co <c 0 + A(Q(T))t. (2.19) 

Combining the estimates in all these steps, we then derive the desired estimates A f(T) < Co 
for some T sufficiently small but independent of e and a. Note that the Taylor sign condition 
and the condition that d>(t, •) is a diffeomorphism can be easily justified due to our estimates 
of time derivatives. Finally we remark that for each fixed a > 0, the Taylor sign condition is 
no longer needed for the inviscid limit problem. This can be seen from (12.71) and (12.81) : even 
g — dfq < 0, one can use the Sobolev interpolation to get the estimates of h for each a > 0. 

We will set the conventions for notation to be used later. N = {0,1, 2 ,... } stands for the 
collection of non-negative integers, N 1+d = {a = ( 0 : 0 , 01 ,..., o^)} is used to emphasize that the 
0—index term is related to temporal derivatives, while is related for spatial derivatives. For 
a € N 1+rf , Z a = <9“° Z^ 1 • • • Z^ d . The Einstein convention of summing over repeated indices will 
be used. Throughout the paper C > 0 will denote a generic constant that does not depend on 
the data, the surface tension coefficient a and the viscosity coefficient e, but can depend on the 
other parameters of the problem, g,n, m > 14 and fi. We refer to such constants as “universal”. 
Such constants are allowed to change from line to line. We will employ the notation a < b to 
mean that a < Cb for a universal constant C > 0. Throughout the paper, the notation A(-, •) 
stands for a continuous increasing function in all its arguments, independent of e and a and 
that may change from line to line, and Ao = A(^). 

The rest of the paper is organized as follows. We collect some analytic tools related to 
Sobolev conormal spaces, the properties of Poisson extension and some geometric estimates 
in Appendixes ei m and [Cl respectively. In Section [3] we study the equations satisfied by 
(Z a v, Z a q, Z a h) and present the estimates of the commutators. Section U is devoted to derive 
the pressure estimates using elliptic regularity in Sobolev conormal spaces, and Section [5] con¬ 
tains the smoothing regularity estimates of h due to viscosity and surface tension. In Section 
[6l the conormal estimates of the solution are derived, and the conormal estimates for normal 
derivatives are given in Section [3 In Section [31 we prove the needed L°° estimates for normal 
derivatives. Finally, the proof of Theorem 12. II is given in Section [9j 
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3. Equations satisfied by (Z a v, Z a q , Z a h) 


3.1. A commutator estimate. In order to perform higher order conormal estimates, one 
needs to compute the equations satisfied by (Z a v, Z a q, Z a h), which requires to commute Z a 
with each term in the equations (11.91) . It is thus useful to establish the following general 
expressions and estimates for commutators to be used often later. 

We will not commute Z a with df directly. For i = 1, 2, 3, set 

Z°dff = df Z a f - dffdf Z a rj + Cf (/), (3.1) 

where the commutator Cf{f) is given for a/0 and j / 3 by 


W)=c?!(/) + c? 2 (/)+c? 3 (/) 


with 


not 

H 1 “ 


nd 

K 2 - 


net _ 

3 — 


z a ,§^,d z f 

OzV 


-d z f 

dj<p 

d z <p 


Z a ,di<p , 


1 


d z ip 


- diipd z f 


Z° 


’ (WJ 


dzV, 


[z a ,d z \f + 


dup 

(9M‘ 


d z f[Z a ,d z \r], 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


for any \a'\ = 1. Note that for i = 1, 2 , ditp = d^ and that for a / 0, Z a d z ip = Z a d z r]. For 
i = 3, similar decomposition for the commutator holds (basically, it suffices to replace ditp by 1 
in the above expressions). Since df and df commute, it holds that 

Z a dff = df {Z a f - dffZ a rj ) + dfdffZ a r] + Cf{f). (3.6) 


It was first observed by Alinhac [ 2 ] that the highest order term of q (which is difficult to control 
for the case a = 0) will be canceled when one uses the good unknown Z a f — df fZ a r), which 
allows one to perform high order energy estimates. 

Since the expressions as f /d z (p will appear often later, we shall first state a general estimate. 
It is assumed that d z ip> and |/i| 2 ,oo < 


Lemma 3.1. For every k E N, it holds that 


dzV 



Proof. Since d z cp = 1 + d z r 7 , so 

f f f d z r] 

d z p 1 l + d z r] 



= f ~ fF{d z rj), 


(3.7) 


where F(x) = x/ (1+x) is smooth and bounded together with all its derivatives on 1+x > cq > 0 
and F( 0) = 0. Consequently, the product estimate (IA.1D implies that 

< ll/llx* + ll/llx* \\F(d z rj)\\ Y , + ||/|| ¥ , \\F(d zV )\\ xk . 

Notice that 

lim>7)ll Y * < A (i ||V,|| t4 ), 

and (1A.1I) implies again, 

ll-^MIIx* < A Q-, ||Vr/|| Y ^ \\d z rj\\ xk . 

Hence the estimate (13.71) follows from these, (IB.71) and (IB.61) . □ 


d z y 


Next lemma deals with the estimates of the commutators Cf(f). 
































10 


YANJIN WANG AND ZHOUPING XIN 


Lemma 3.2. For 1 < |ce| < m, i = 1, 2, 3, it holds that 

1 


\\C? (/)|| < A \h\ Y ^ +1 + \\d z f\\ Y ^j {\\d z f\\ xm ^ + \h\ xm _ h i ) . (3.8) 

Proof. We only present the proof for i = 1, 2, and the last case is similar and slightly easier. 
First, for Cf t , it follows from the commutator estimate (1A.4I) that 


\r a II < 

Wl lb 


dup 


d z ip) 

Consequently, (13.71) yields that 


_ 2 \\Zd z f\\ xrn - 2 + \\Zd z f\\ Y ^2 


z 


djip 

d z ip 


fell < A ^ —, ||V<^|| Y >™ + |/i| ¥ mm +1 + \\d z f\\ Y ™.^J ^||^/|| Xm ^i + |/i| x9 n_ 1; i + llc^llx™- 1 
It then follows from (11.51) . (IB.71) and (IB.61) that 

fell < A \h\ Y v +1 + \\d z f\\ Y ™) (||^/|| xm _x + \h\ xm _^ ) . (3.9) 

Next, for the first term in C“ 2 , one can use similar arguments: (IA.4I) and (11.51) yield that 


d z f 


Z a ,dnp, 


1 


d z ip. 


< 


/II L°° ( \\ZdiP\\xm-2 


Zd z ri 


{dzPf 


+ ll^j^H 


Zd z r] 




and hence by using (13.71) . (IB. 61) and (IB.71) . one can show that 

1 


dzf 


Z a ,di<p, 


d z tp 


< A \h\ Y y+i + \\d z f\\ LOO '] 


By (IA.2I) instead of (IA.4I) . the same estimate holds for the second term in C“ 2 . Hence, 


|C“ 2 || < A ( —, \h\ Y ™ +1 + ||9~/|| L o 

, c o 


00 / I * ^ I w-m— 1 


1 . 
2 


It remains to estimate C” 3 . Notice that 


[Z a ,d z ]f= ]T cpd Z (zPf) 


'\<m -1 


for some harmless smooth bounded functions eg. This yields, by using again (IB. 61) . 

fell < A (j- } \\di<p\\ Loo + \\d z f\\ L J) {.\\d z f\\ X m-l + IIMxm-O 


< A \h\ loo + \\d z f\\ Loa '] (||3 z /|| xm _i + • 


Consequently, the estimate (13.81) follows by collecting (13.91) . (I3.10|) and (13.121) . 


(3.10) 

(3.11) 

(3.12) 

□ 


3.2. Interior equations. We shall now derive the equations in the domain 11 satisfied by the 
good unknowns V a = Z a v — dfv Z a rj and Q a = Z a q — dfqZ a rj. 


Lemma 3.3. For 1 < |a| < m, it holds that 

dfV a + v • + V^Q" - 2 eV v • S v V a 

= (dfvV‘ fi v)Z a ri-C a {T)-C a {q)+£V a (S‘ p v) + • (S a (v)), 

V‘ p -V a = -C a (d), 

where the commutators C a (q), C a (d), 8 a {v) andC a (T ) satisfy the estimates: 


\\C a (q)\\ < A (1 |fc| Y * + i + ||d 2 g|| Y ™ J (||^9||xm-i + , 

||C"(d)|| < A Q-, \h\ Y rn +1 + ||0 a u|| ¥ ^ (||9 z u]| xm _i + |/i| xm _!,i ) , 


(3.13) 

(3.14) 

(3.15) 

(3.16) 
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»£»ll<A(i |ft| Y » + ‘+ (||d ; 




\rrn — 1, tt J 9 


(3.17) 


||C“(T)|| < A \h\ Y rn +1 + ||w|| y¥ + ||a z u|| Y ™^ (||u|| x ™-i + IIMx^-i + > 

and T> a (S v v) is given by 


(3.18) 


V a (S*v) ij = 2Cf(S'Pv) ij . 

It should be noted that the commutator 1^(S^v) will be estimated later by using the inte¬ 
gration by parts. 

Proof. First, the equations (I3.13D — (|3.14[) follows from applying Z a to the equations (11.91) . In¬ 
deed, (13.61) implies that 


Z a V^q = + dfV v qZ Q rj + C a (q), 

where C a (q) = (Cf (g),Cf (q),C%(q))\ and 

Z a V 1 ^ • v = ■ V a + dfV* ■ vZ a r] + C a (d), 

where C a (d) = Yll=iC?(vf). 

Next, note that 

df + V ■ = dt + Vy ■ VyV + V Z 3 Z , 

where V z is defined by 

V z = — — v z with v z = v ■ N — dt<p = v • N — d\rj. 
d z ip 


(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 


By using (13.211) . one can thus get that 

(df + vV f )v = (d t + vyV y + v z d z ) Z a v + (u • Z Q N - d t Z a r])dfv 
- dfZ a rj(v • N - d t r])dfv + C a (T) 

= (df + v ■ V^)Z a v - dfv(df + v ■ V^)Z a rj + C a (T ) 

= (df + v-V^)V a + df (df + v ■ V?) vZ a r 1 - dfv ■ V^vZ a v + C a (T), 
where the commutator C a (T ) is defined by 


C a (T) = [z a , Vy\ dyV + [Z a , V z , 9 z v] + 


z a ,v z , 


d z ip 


d z v + — [Z a , v} • N d z v, (3.24) 

C'zP 


+ v z d z v 


Z c 


1 


Z a 'd z r, + V z [Z a , d z \ v + [Z a , d z ] V 


’ (d z p) 2 \ ~ ’ (d z vf 

for any \ot'\ = 1 . 

It remains to compute eZ a /SPv = 2 eZ a V ,f ■ (S^v). Note that 

2 • (S^v) = 2V^ • (Z a S^v) - 2 (df S tp v)V fi (Z a r ] ) + V a (S*v) 
with T> a (S tp v) i = 2Cf(S‘ f v)ij, and 

2 Z a (S*v) = 2 S* (Z a v) - dfv < 8 > V v Z a r) - \/ v Z a r] <g> dfv + £ a ( v ) (3.25) 

with (, £ a (v ))^. = Cf(vj) + Cf(v t ). Hence one deduces that 

sZ a A^v = 2 e • S v (Z a v) - 2eV* • (dfv (8) - V v, Z a 7 1 ® dfv) (3.26) 

- 2e(df S‘ f v)V v (Z a r ] ) + eV a (S lp v) + eV^ • (£ a (v)) 

= 2 e • S v V a - 2 e(dfV p • S v v)Z a r] + eV a (S^v) + eV^ • (£ a (v)). 


Consequently, these and (|1.91) imply (13.131) (13.141) . 
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Now we estimate these commutators. First, thanks to Lemma [3. 2 1 the estimates (I3.15I) (I3.17D 
hold. To estimate the commutator C a (T) defined by (13.241) . one needs to bound v z and V z . It 
follows from (IB.71) that 


Klly? + - A > IMI Y ? + ll Vr ?ll Y i + II^IIy^) - A IMI Y ¥ + • 

(3.27) 

And (1A.1D . (IB.61) . and (IB.71) yield that 

\\Zv z \\xrn-2 < \\Zdtr]\\ X ru-2 + ^1 + ||Vl?|| y ,m- 2 ^ |M| X m-l + IMIy 7 "^ 2 ll^llx^- 1 (3.28) 
< A (||v|| y ™^2 + |h| ¥ -) (|M| X m-l + • 

With (I3.28P and (13.271) in hand, by using (IA.1D . (13.71) . (IB.61) and (IB.71) . one can obtain that 


\ZV~ 


< 


f 1 \ 


1 

Z — v z 

+ 

Zv z 

\dzpj 

X m-2 

d z <P 


(3.29) 



Consequently, one may use (IA.2I) . flA.4p . (|3. 1 111 and (13.71) combined with (13.271) (I3.29P and also 
again (IB.6P , (IB.71) , similarly as in the proof of Lemma 13.21 to conclude the estimate (13.181) . □ 


3.3. Boundary conditions. We shall now also compute the boundary conditions satisfied by 
(Z a v, Z a q, Z a h) when 03 = 0 (for 03 7 ^ 0, Z a v = 0 on the boundary). As a preliminary, one 
has one has the following. 

Lemma 3.4. For k € N: 

|Vu| X fc < A , |/i |^| +1 + ||Vv|| y fc^ (Mxm + |h| X fc,i) (3.30) 

and 

|Vu| xM < A \h\ Y k +2 + ||Vu|| y ! +1 ^ (|u| xfclS+ i + |h| X fc, s+ i) for s = (3.31) 

Proof. Note that it suffices to prove the estimates for d z v. Since ■ v = 0, thus 

d z tp (<9iui + <9 2 u 2 ) + d z v • N = 0. (3.32) 

Then for s = — ||, (IA.12D implies that 

\d z v ■ n| X fc, s < A , |Vry| Y fe +1 + ||V 3 / u|| y fe +1 ^ (MxM+i + l^hlx fc > s ) (3.33) 

~ ^ (co ’ |/l| Yl+ 2 + H V 2/ U II Y §+1^ (HxM+i + WxM+i) ■ 

where the second inequality follows from Lemma IB. 21 and the trace estimate (IA. 6 D . 

To bound T\d z v, we shall use the boundary conditions in (11.91) which yield 

n (S^vn - n X v) = 0. (3.34) 

To compute II(S' ¥, un), one can use the local basis (d y i , d y 2 , d y 3 ) in Q t induced by (11.41) . The 
induced riemannian metric is given by g z j = d y i ■ d y j , whose inverse denoted by g lJ . It follows 
from the definition and (11.81) that (d y iu)(t,&(t, ■)) = dtv. Hence, 

2Sun = n • Vu + Vmj= d n u + g ZJ d y jU ■ nd y i = d n u + g l] djv ■ n d y i. 

Note also that 


(3.35) 
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one then gets from (13.3411 that 

r\ 

Tld z v = — (di^plldiv + d‘ 2 iplid 2 v — g^djV ■ nil d y i — n-xfilv) . (3.37) 

Hence, using (IA.12I1 again shows that 

|n<9 2 u| x /c, s < A + ||Vn|| Y | +1 ^ (M x m+i + |h| X fc, a+ i + \d z v ■ n| x ^ jS ) (3.38) 

< A \h\ Y § +2 + ||Vn|| Y | +1 ^ (M xfc , s+ i + |/i| X fc, s+ i). 

where in the second inequality (I3.33H has been used. 

Consequently, we conclude the estimate (I3.31H by combining (13.3311 and (|3.38l) . And (13.3011 
follows similarly by using (lA.lH instead of (1A.12H . □ 

We now study the dynamic boundary condition on {z = 0} and the Navier slip boundary 
condition on {z = —b}. 

Lemma 3.5. For 1 <\a\ <m such that C 13 = 0, it holds that on {z = 0} 

2e5 ¥, H“N- (Z a q-gZ a h + crZ a H)N (3.39) 

= - 2 e^unZ Q N - 2eZ a hdf(S‘ fi v)N + eC a (B £ ), 
where the commutator C a (B £ ) satisfies the estimate: 

|C a (£b)lo < A ^ —, |/i| Y ™+i + ||Vu|| Y m^ (M X m-i,i + |h| Xm _i,i). (3.40) 

Similarly, on {z = —b} one has that 

V? = 0 , (S^V a e 3 )i = KVfi - S a (v) i3 , i = l,2, (3.41) 

and 

\£ a (v)\ 0 < A , \h\ Y in +1 + \\Vv\\ Y ^j (M X m-i,i + |/i| xm _i,i). (3.42) 

Proof. Applying Z a to the dynamic boundary condition and using (13.251) . one gets 

e (2( Z a v) - dfv (8) V v Z a ip - V v Z a v <g> dfv + e£ a {y)) N - (Z a q - gZ a h + aZ a H ) N 
= - (2 eS*v — (q — gh + oH)I) Z“N - [Z a , 2 eS v v - (q - gh + aH)I, N] . 

= — 2e (S v v - S^v n • nJ) Z“N - 2£ \Z a , S^v - S^vn • n I, N] (3.43) 

= —2eS lf vIlZ a N - 2e [Z a , S*v n, N]. 

This yields (|3.39l) with the commutator C a (B e ) defined by 

C a (B £ ) = -£“(u)N - 2 [Z a , S^vU, N] . (3.44) 

Similarly, applying Z a to the Navier slip boundary condition and noting that V a = Z a v, one 
shows (13.4111 . 

Now, it follows from (1A.4I) . (IB.7H and Lemma l3~41 that 

|[z“,5 ¥ ’un,N]| 0 < \z(s^vU)\ xrn . 2 |ZN| Y m_a + \z(s*vn)\ Y rn ^2 |ZN| xm _ 2 ( 3 . 45 ) 

< A l|Vr/|| ¥ - + ||Vn|| Y -^ (|Vu| X m-i + |/»| X m-i,i) 

< a (^’ + II Vv IIy^) (Mx—m + Wx-m) ■ 

On the other hand, following the proof of Lemma 13.21 and using again Lemma 13.41 one has 
\£ a (v)\ 0 < A Q-, llVr/H^ + ||Vn|| Yl ) (iVi^-, + \Vrj\ X m-i) 

< A ^ —, \h\ Y rn +1 + ||VU|| Y -) (|u| Xm _l,l + |/l| X m-l.l) • 


(3.46) 
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Consequently, the estimates (13.4211 and (13.401) follows. □ 

Finally, we study the kinematic boundary condition on {z = 0}. 

Lemma 3.6. For 1 <\a\ <m such that ot% = 0, it holds that on {z = 0} 


d t Z a h + v y ■ V y Z a h — V a • N = -dfv • N Z a h + C a {h) (3.47) 

where the commutator C a (h) satisfies 

\C (h)|g < A , \h\^rn +1 + ||u|| Y m^ d^lx" 1-1 - 1 + Mx™- 1 ) • (3.48) 

Moreover, 

C a {h)= Y C^Z a ~ a 'v y ■ S7 y Z a 'h + C Q (h), (3.49) 

l«'l=i 

andC a (h) satisfies the estimate: 


C a [h) < A —, |/i| m +2 + |v| m +1 


co 


i-l,2 + C 


Proof. Applying Z a to the kinematic boundary condition yields 

d t Z a h + v y ■ V y Z a h — Z a v ■ N = C a (h) 

where 


(3.50) 

(3.51) 


C a (h) = -[Z a ,v y ,Vyh\. (3.52) 

This yields (13.471) . and the estimate (13.481) follows from (IA.4D . One may further single out the 
highest order derivative terms according to (13.491) . where C a (h ) is defined by 

c°(h) = - £ CfaZPyy • VyZ^h. (3.53) 

/ 3 + 7 = a 

13 ^ 0 , h \>2 

And the estimate (13.501) follows by using (IA.1D . □ 


4. Pressure estimates 

In view of the equation (13.131) . one needs to estimate the pressure q. The first equation in 
(11.91) implies that 

A = -V* ■ (v ■ V^v) in fl. (4.1) 

Moreover, the dynamic boundary condition gives 

q = 2 eS^v n • n + gh — crH on {z = 0}. (4.2) 

Projecting the first equation in (11.91) along N onto {z = 0} and {z = —b} yields 

W^q • N = -d t v • N -(vy V y )v ■ N + eA^v ■ N on {z = 0} (4.3) 

and 

V^q • N = eA^v • N on {z = -b}. (4.4) 

Here in (14.41) one has used the fact that N = e^ and V 3 = 0 on {z = — 6}. 

Note that to solve the pressure, one has two choices of boundary conditions on {z = 0}, 
he., <m and (14.31) . Without surface tension, one can use the elliptic problem (14.11) . (14.21) and 
(14.41) to establish the regularity estimates for q. The subtlety lies in that the energy dissipation 
estimates of (USD in the case without surface tension provide the needed estimates for those 
boundary terms. When there is surface tension, however, the energy dissipation estimates do 
not provide enough estimates for the boundary term — crH (which is of one half regularity less). 
This would suggest that the elliptic problem (14.11) . ()4.2|) and (14.41) is not the right choice for 
estimating the pressure q in the case with surface tension. Our way to get around this difficulty 
is to use instead the elliptic problem (ED, (USD and (ED- It is then noticed that this approach 
forces one to estimate the time derivatives of v, that is, one needs to perform energy estimates 
for the time derivatives of the solution. However, there is an essential difficulty arising: when 
doing energy estimates with time derivatives up to m order, we can only obtain the estimates 
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of time derivatives of q up to m — 1 order due to the presence of dfi) in (14.31) . Thus the energy 
estimates cannot be closed since it seems that m order time derivative of q is involved. We will 
explain this and our way to overcome it in more details in Section 16.31 
It follows from the definition of df that 


X* ■ v = -J—X • (Pv), X v f = P = 


d z ip 

And then ISP can be expressed as 


d z <p 


dzV 

0 


0 0 

d z tp 0 


—d\ ip -d 2 tp 1 


A*7 = • (V*7) = — V • (EXf) 


(4.5) 


(4.6) 


with the matrix E defined by 


E = 


d z ip 


d z (p 0 —d\ <p 

PP* = | 0 d z ip -d 2 cp 

-d W -d2<p i+(ai i 2 : (82y)2 


Note that E is symmetric positive and that if ||Vj,^||ioo < A and d z ip > co > 0 then there 
exists 5(co) > 0 such that 

EX ■ X >6\X\ 2 , MX € 

Moreover, 

1 


|£|| ¥ * — A | ^ , j/llyfe + l ) • 


(4.7) 

(4.8) 


One can write 


if = Id + if, E = 


d z r] 0 —d\ rj 

0 d z r] ~d 2 r] 

-dm -d 2V ‘.omr+iihf-s.v 


where 


E 


1 


<A l-,\h\ 




(4.9) 


To' 1 'y* +1 . 

Here t k and Y fc are referred to the usual spatial-time Sobolev spaces as defined similarly as 

(HU). 


Since N = P*e 3 on {z = 0, —b}, the equations (14.11) (14.41) can be rewritten as 
- V • (EXq) = F := d z ipX^v • X*v in Q, 

q = G 1 := 2 eS v v n • n + gh — aH on {z = 0}, 

EXq ■ e-s = G 2 := —dtv ■ N — (v y ■ X y )v ■ N + e/S^v • N on {z = 0}, 

EXq ■ e-s = G 3 := eA^v • N on {z = —b}. 

We shall now prove the estimates for the pressure q. 

Proposition 4.1. The following estimates hold: 

Ikllx* HVgllx* P ll^z^llx' 1 - 1 

<A(—, \h\ k ,, + \h\ fc +7 + Ilf II fc +3 + HVf II k, 2 + ||f || fc+7 + ||Vf|| k ±5 

— ^co 'y ^+ 3 1 'x“2- 11 "yt 11 m y ^+ 2 11 "xT- 11 "x-f 

x (|ii| xfc+li _ 1 +a\h\ X k ,2 + ||u|| X fc+i + ||Vf|| xfe + £ |f| xfc ,3 + e N x fc,i|) » fork> 3; 

Ikllyfc + II v 9|Iy* + ll^^llYfc- 1 (4-15) 

< A l^l ¥ fc +4 + \h\ X k+4 + |HI Y i+5 + ||Vw|| y * +3 + ||w|| X fc+4 + ||Vu|| xfc+3 ) , for k > 1. 


(4.10) 

(4.11) 

(4.12) 

(4.13) 


(4.14) 
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Proof. Multiplying the equation (14.101) by q and then integrating by parts over 17, using (14.121) 
and (|4.13l) . one obtains 

(EVq, Vq) n = (F. g) n + (G 2 , q) z=Q - (G 3 , q) z= _ b . (4.16) 

Here (•, -)q , (-, -) z=0 and (■, -) z __ 6 denote the L 2 inner products on 17, {z = 0} and {z = —b}, 
respectively. It follows from the trace estimate |g| i ^ Iklljfi and Cauchy’s inequality that 

||Vg|| 2 < \ (||F|| 2 + IG 1 ^ + |G 2 | 2 _i + |G 3 |\) + r, \\q\\ 2 Hl (4.17) 

for any rj > 0. By the Poincare inequality (1A.8I) and (14.111) . one can get by taking rj sufficiently 
small that 

h\\m — Ao (||-F|| + |G fl | 0 + |G 2 |_i + |G" 3 1 i) • (4-18) 

Note that if one uses solely the problem (14.11) . (14.21) and (14.41) to estimate p, then one needs 

G 1 i. This half less regularity requirement enables us to control the surface tension term by 
2 

the energy dissipation estimates. 

Next, applying Z a with |a| = k to the equation (|4. lOj) and using (14.121) (14. 13|) lead to 

(Z a (E ■ Vg), VZ“g) n = ( Z a F , Z a q) Q + (Z“G 2 , Z a q) z=Q - (Z a G 3 Z a q) z= _ b . (4.19) 

Then as for (14.181) . one can derive after using (14.111) that 

Ikllx* + HVgllx* < A 0 (ll^llx* + l^lx* + I^Vi + |G 3 | x fc,_i (4.20) 

+ \\E-[Z a ,V]q\\ + \\[Z a ,E]-\7q\\y 

We then estimate the commutators in the right hand side of (14.201) . First, (13.111) implies that 

II E • [Z Q , V] g|| < \\E\\ Loo ||Vg|| x *-i < W ¥ i || Vg|| xfc -i • (4.21) 

However, the other communtator needs more attentions according to 1 < k < 3 or k > 4. 

Indeed, for 1 < k < 3, direct estimates by controlling the Vg terms in L 2 and E terms in L°° 

yield 


II [Z Q ,E] • Vg|| < \\E\\ fk || Vgllxfc-i < \h\ Yk+1 ||Vg|| x *-i • (4.22) 

Plugging (I4.21|) - (I4.22|) into (14.2UD . by an induction argument and (|4.181) . one can deduce that 


+ ll^^llx^ — ’ I|£|Iy3 

Co 


+ IG^ 1 1 u + |G 2 | ,_1 + |G 3 


F-h + 


< A 


1 

co' 


+ \G 3 L k + |G 2 


+ G 3 


Since the equation (14.101) gives 


d zz q = {F-d z | E 3d d j<l ] - di 


i<3. j 


in a similar way and by (|4.23|) . one can also obtain 

1 


l|3**g||x*-i < A -ll^ll 
To 


1 


<A —,|*h 

To 


-i + 

+ \G 3 L k + | G 2 1 


+ G 3 


It then follows from (14.23|) and (|4.25|) that for 1 < k < 3, 

+ HVgllx* + ||^g|| X fc-i 


< A 


(4.23) 


(4.24) 


(4.25) 


(4.26) 


\ - , \h\ Y 4 + | h\ k+ 3 + \h\ fc+3 1 ~f" 11 F 11 fc+3 T I G^ | k+ 3 + I G 2 | fc+3 1 + I G 3 | fc+3 1 ] 

^Co WT 1 'xT-’ 5 " 11 XT 1 X 2 I 'XT-I 1 'XT "2 y 
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x ( Wx fc >2 + + \ Gl \ x * + + I g3 | x *.-0 • 

We now claim that (|4.26f) holds for all k > 1. This will be proved by induction. Assume that 
k > 4 and that (14.261) holds for k — 1. The commutator estimate (IA.2D yields that 

||[2T,£] • Vg|| < ||Vg|| ¥ ^i + WZE^ ||Vg||x*-i (4.27) 

< A N y *+3 + ||Vg|| YV ) (|/i| xfc ,i + ||Vg|| X fc-i) . 

Plugging (14.271) and (14.211) into (14.201) and using (14.241) again lead to 

IMIx* A ll v 9llx* 4“ 11 ^zz Of11 x^ 1 —A , |h| Y fc+3 + HVgl^fc^i ^ ^ |^l x fc,i 4" ll^ ( lllx fe ~ 1 (4.28) 

+ l|A|lx fc + | g1 |x* + l G2 |x fe ’-2 + l G3 lx‘.-i) ‘ 

To remove the dependence of A on ||Vg|| k- 1 , one applies the anisotropic Sobolev embedding 
estimate (IA.5D and the induction assumption to obtain 


l|Vg|| Y ^i < ||<9 2 Vg|| x fc+i ||V g || x ^ (4.29) 



Here one has used the fact that AiH < k — 1 since k > 4, which allows one to use the induction 
assumption. Plugging the estimate (|4.29l) into (|4.28l) and using the induction assumption to 
estimate ||V(?||xfc-i, one thus concludes (J4.26P for all k > 4. 

We now estimate the right hand side of (14.261) for k > 3. It follows by the product estimate 
(IA.1I) that 

ll-F’llx* = V^llxk (4.30) 

SA^.IM^ + IIWIIyi) (l*U + I|W|| X ,). 

Similarly, since k > 3, 

ll^ll x ^ ^ A , |/i| y a^3 +1 + ||Vu|| y ^) i + ||Vu|| x ^) (4.31) 

A ( — , | /i| k +1 + \h\ k +2 4- ||Vu|| k T ||Vu|| fc+3 ) • 

— ^co 'y? + 1 1 'x^+ 2 11 "y^ 11 "x-2-y 

By (IA.1D and Lemma 1X41 

\G\ k = |2eS*W n + gh-aH\ xk (4.32) 

< A , \h\ Y k +2 + l|Vu|| Y fc^ (e |Vu| xi! + |h| X fc + cr|^|x fc . 2 ) 

< A , \h\ Y k +2 + ||Vu|| y ,) (e Mxm 4- \h\ X k + & H x m) ■ 

Similarly, since k > 3, the trace estimate implies that 
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- A (^’ |/l| Y l + 2 + + H Vw M + IHI x -§+ 3 + l|Vv|l x§+ 2 ) ' 

To estimate the most delicate term G 2 , we start with eA^v ■ N. Note that 
A^v • N = 2 (V v ■ S*v) • N = 2V^ • {S^v N) - 2 S v v : V^N 


and 

• (S^uN) = di (S^v N)j + d 2 (5 ¥, uN) 2 + -^d 2 (S^uN) • N. 

OzV 

Hence by the estimate (I A. 121) with s = — 4, 4, one can get 


But 


\A*v N| 


< 


d z p 


d 2 (S^uN) • N 


x + + |S*t, : V^Lj 


< A \h\ Y * +3 + ||Vu|| y|+1 + \\d z (S*v N) ■ N|| y * +1 ) 
x (\d z (S^vlSS) ■ N| x j. _i + |Vu| xi i + |^l x fc,|) ■ 


<9 2 (5^uN) • N = d z {S^v N • N) - SOiN • d 2 N, 
and recalling the matrices P and E, 

d z {S^v N • N) = d z (V^uN • N) = d z (^—P*VvP*e 3 ■ P*e 3 J 

\o z p J 

= V^(9 2 u)N • N - d jVi d z (E 3l P 3j ), 
one can then deduce from (14.361) that 

I • N| x „_j < A (±, |A| y ,„ + II Vv|| ¥ , +1 + l|W(A»)N • N|| y4+1 ) 

x(|W(9 2 »)N.N| xl ,_ J+ |VH xl , 4 +|/.| xfcJ ). 

Note further that 


(4.34) 

(4.35) 

(4.36) 


(4.37) 

(4.38) 

(4.39) 


V v (<9 2 w)N • N = V* 3 {d z v • N) • N — (N • V v ) N • d z v, (4.40) 

and one can compute by using ■ v = 0 that 

{d z v • N) • N = -V v ( d z (p(d lVl + d 2 v 2 )) ■ N. (4.41) 

It follows from these, (14.391) and Lemma 13.41 that 

IA'VNI^j < a( 1 IM y4+ , + I|V„|| y4+1 ) (|Vi.| xfc , +|fc| xt .) (4.42) 

£ A ( x - + l|w|l vi+>) 0" A + ■ 

Next, one easily has 

IK • V> • N| xfci _i < A Q-, \h\ Y , +2 + ||v|| y4+2 ) (|w| xfcli + l^l xfc ,i) • (4-43) 

Finally, we estimate the remaining time derivative term dtv ■ N. One first has 

\d t v • Nl^.^i < A Q-, \Vh\ Y k-i +1 + ||5 t u|| ¥ fc-i +1 j (4.44) 


< A (1, |h| Y ^3 + IMI Y l+a) (\v\ xk i + • 

It then suffices to estimate \d^(dtv ■ N)|_i. However, this will lead to some difficulties since k 

2 

can be m — 1, and energy estimates yield only d™v € L 2 (H) which cannot ensure the control of 
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the H~z{z = 0} norm of d™v. The key observation is that • N is indeed in H~z ({z = 0}). 
The way of achieving this is to use the Alinhac good unknown, returning back to k, V k+1 = 
d k+l v — dfvd k+1 rj. Indeed, since • V k+l = — C k+1 (d ) with C k+l {d) dehned as in (13.1411 with 
Z a = d k+1 , by Lemma TC.31 and using the estimate (|3.16l) and (IB.pp . one gets that 


V k+1 • N 


< A 
1 ~ 

' 2 

< A 


s)( 


yk+l 

> II Vn || L 


+ 


. yk+l 


(4.45) 


9f +1 „ 


+ 


d k+1 rj 


+ 


C k+ \d ) 


<a(L,|/.| yS ji +i + ||v»ii y4 j 1 ) (||af +1 . 


|ft|. 


: fc + i ’-2 


This in turn implies 


d$ +1 v • N 


< 

1 ~ 
" 2 


< A 


v k+1 • N 

1 


+ 


N • dfvdf+'r) 


(4.46) 


, \h\ fc + 3 + ||Vv|| fc+1 

c 0 'f-r 11 ¥ 2 ~ 


a fc+1 ? 


Hence, it follows from (I4.44[h (14.46p and the trace estimate 
I dtv • N| i < \d t v • N| fc _j i + d k (d t v • N) i 


+ I^Lfe+i.-i + 


< I • 11 < IN 


m 


that 


(4.47) 


< A W y *l+ 3 + Ill’ll y 4+4 + HVcH^fc+i J (JM| x m-i + 11Vi;11 X fe + |/i| xfc+1 _ 

Therefore, we conclude from (14.421) . (I4.43p and (14.471) that 


\G 2 


'.-h ^ a Q)’ i^i ¥ ! + 3+ m y h¥ +ii vvii Yi+a ) 


(4.48) 


x (ll' u llx fc + 1 + HVullx* + |/l| x fc +1) _i + £ |^l x fc,3 + £ |^l x fc,3 
Similarly, since k > 3, due to the trace estimates, one can get 


G 2 \ fc+3 1 < A ( , \h\ fc+3 , o + ||u|| fc+3 . o + || Vn|| fc+3 


Co ’Y 


i+2 


(4.49) 


X ( 11111| fc+3 I , T || fc+3 + \h\ fc+3 . . 1 — f- £T I 'll | fc+3 3 + £ | h\ fc+3 3 ) 

V XT + 1 XT IT H -2 1 'XT’2 1 'XT'd 


— ^ N y §+3 + l^l X T + ll Vw ll Y §+2 + IMI X *+1 + l|V?;|| x fc+B^ . 

Consequently, plugging the estimates (14.301) (14.31 1) . (|4. 321) - (14.331) and (14.481) (14.491) . along 
with doing the the same estimates for G 3 as for G 2 , into (14.261) . we can obtain that for k > 3, 

Nix* + l|V?|| x * + II^NIxfc-i (4.50) 

— ^ (eg 1 ^Y2+ 3 ^^X^T H*Hyt II^' W IIy2+ 2 II^IIxt^ 

X (|/i| x fc +1 ,_ i + cr\h\ X k ,2 + ||c|| X fc+i + ||Vi;|| X fc + e M x fc,3 + e l^l x fc,§) • 

This proves the estimate (14.141) . 

To prove ()4. 151) . one can use the anisotropic Sobolev embedding estimate (1A.5D and the trace 
estimates to have that for k > 1, by (14.141) . 

IMIy* + l|Vg|| Y * ~ l|Vg||x fc+2 + ll^ 22 (?|lx fc + 1 ( 4 -51) 

— f —i \h\k +4 4- \h\ fc+9 + ||n|| fc+s + ||Vn|| fc , 3 + ||d| fc+9 + ||Vn|| k+1 

\Co 'Y^+ 4 1 X 2 11 YT 11 "YI +J " X 2 11 X 2 


/■fc+3,— 4 + a 


>2 + IMIx fc + 3 + l|Vc|| X fc +2 4- £ |c| xfc+2 ,3 +e|h| xfc+2i 3^ 


x 
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< A 


co 


;\h\ 


Y2 


+4 


+ 


+ ||V^|| ¥ fc +3 + ||u|| X fc+4 + ||Vp|| xfc+ 3^ 


Note that one has used the fact that k + 2 > 3 so that (|4.14l) can be used with k + 2. This, 
together with (|4.24[) again, proves (|4.15l) . □ 

5. Smoothing estimates of h 

We first show the smoothing regularity estimates of h coming from viscosity. 

Proposition 5.1. For every m € N, e € (0,1), it holds that 


£ I^)I x ™-i,3 <e\h(0)\^ n _ lt i+ J A 


k + 2 + ( £ ML-i,3 +£\h\ 2 


m-l,| 


and 


(5.1) 


£\h(t)\ 2 i < s\h{0)\ 2 i + / A 


+2 £ |u| 




(5.2) 


Proof. We prove only the estimate (15.21) , and the estimate (15.11) follows in the same way. Apply 
Z a with a € N 1+2 , | a\ < m, to the kinematic boundary condition to get that on {z = 0}: 


d t Z a h + v ■ X7 y Z a h - Z a v 3 + [Z a , v y ] ■ V y h = 0. 

Then applying further A 2 , the tangential Fourier multiplier, to (15.31) gives 

d t AFZ Q h + v -V y hhZ a h - A^ Z a v 3 + [A3,uJ • W y Z a h + A^ ([Z a ,v y ] • \7 y h) = 0. 
A standard energy estimate on the equation (15.41) yields 


d 


— \Z a h\\ < ||V„u|| roo \hT 1 
dt 1 b ~ 11 y UL 1 


+ 1 Iwl xm ,i + 


A 2 ,u, 


• VyZ a h 


+ \[Z a ,Vy\-Vyh\l |hl 1 


Due to the commutator estimate (1A.10I) . one has 


A 2 , v y 


■ VyZ°h 


< ||V y u|| i<x , \h\^ 


1 + HV^IIj 


1 . 
>2 


And the estimate 1 A. 121) leads to 


| [Z a ,Vy\ ■ Vyh\l < Yl \Z fi ~ fi, Z^Vy ' V yZ^ h 


I0I + |7| = “ 


< 


m m 

1 Hy’it 


+2 I^L m,l + \h\ Y ^+2 


(5.3) 

(5.4) 

(5.5) 


(5.6) 


(5.7) 


Hence, plugging the estimates (15.61) and (15.71) into (15.51) and summing over |a| < m, by Cauchy’s 
inequality, one can deduce that 


ill h\i,<(l + \h\ 


+2 


- \h \ 2 , 

dt x m ’7 

Integrating the inequality (J5.8D directly in time yields (15.21) . 


YT+ 2 II l ’llY"T 


v\ 2 1 + \h\ 2 1 

1 X m ’2 1 'x m ’2 


(5.8) 

□ 


Next, we show the smoothing estimates of h due to surface tension. 
Proposition 5.2. For every m € N, e, a € (0,1), it holds that 


WL-1.& <a 


1 

co : 


K 


m+3 


+ l|Vu|l 


m -\-1 


2 1 1 + \h \ 2 


m—1, ? 


+£ 2 \v \ 2 , 3 +E 2 \h\ 2 . 3) 


(5.9) 
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Proof. Apply Z a with a € N 1+2 , |a| < m, «o < m — 1, to the dynamic boundary condition to 
have 

- aZ a H = Z a q - 2eZ a (S v v n • n) - gZ a h. (5.10) 

Note that 


Z a H = V„ • Z 



Vyk 


x/i + lV^I 2 


(5.11) 


= V„ 


( V y Z a h 


| = + \7 v hZ a I — 

V VI + IV^I 2 v Ul + I Vyh\‘ 


+ 


z a ,v y h, 




and 




v'l + iv^l 5 


_ gOi — cx. 


1 


\/ y h ■ X7 y Z a h 


yr T\%WJ 

a—a' 


= -z a ~ a 


a-a' ( ^yh ' ^yZ a h 


Wl + IVy/ll 1 


(5.12) 


yi+W 


z c 


Vyk 


\n+W7A‘ 


VyZ° h 


for any \a!\ = 1. Hence, 

Z a H = aV y 

where 

C(H“)= - 


I V y Z a h __ V y fo • v ^ + 

W1 + |V^I 2 Vl + |V y /»|2 3 W 


-1 

N] 

P 

1 

P^ 

< 

• V.„Z a 'hV,,h + 

Z a 1 X7..h 

[ xA + iv^rJ 


[ Vi + W 


(5.13) 


(5.14) 


for any \a'\ = 1. It follows that 

' ■ Ufo “ C(S “ )+z “" - 2£Z " (S ^ ’- n) - 9Z “ ,! ' 


-crV,, 


Then apply further A 2 to (15.151) to get 

( V y hzZ a h V y h-\7yhhz a h 
Wl+IVP ^l + \Vyh[ 2 
1 


- cVi 


V„/i 


■3 v 2/ 


(5.15) 

(5.16) 


= dV, ; 


A 2 . 


- 1 


\7 y Z a h- At. 


^-V,,h- 




’ Vi + iv^l 2 

+ At (crVy • C(B“) + Z a q - 2eZ a {S^v n • n) - gZ a h ). 

It follows from a standard energy estimate for this elliptic equation that 

\S7 yhk Z a k\ 2 \Vyh-VyhhZ a h\ 2 ^ 


3 v v 


VyZ a h 


a 


2=0 \ \A + \^yh\ 2 y/TTWyW 

1 


dy 


(5.17) 


= —o~ 


A 2 


Vi + iv^l 2 


- 1 


V y Z a h~ At. 


V?//l rVjjh- 


x/i + lV^I 2 


■3 v y' 


VyZ a h,\/yA2Z a h 


- a 2 [A^C{B ( f) 1 VyA^Z a h^j + [hh {Z°q - 2eZ Q (^u n • n) - 5 Z“/i), ahhz a h^ . 
Since for any vector a£R 2 , 


|V y /i • a| : 


V^+WW y/1 + \V y hf > vH + IV^I 2 


lal 2 
t 3 l a l 1 


(5.18) 
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so Cauchy’s inequality, together with (1A.11I) and (1A.12D . since |«i| + |CK 2 1 > 1, give 


iVy/iL.! 5 < Aq |g| 1 + e \S^v n • n|; 


+ 


m- 1,7 


(5.19) 


2 


A 2 , 1 = - 1 

V v Z a h 

V 1 + \^y h \ 2 





A*,- 


Vyh 



2 \ 

VyZ a h 

j 


^ A { —, ^Y 12 ^ H Vu ll Y ^ 


Vt+WyW" 

■) (j^X—^ + 

+a 2 |^ m _ 1 ,3+e 2 |^n-n|2 m _ li i). 

Using the similar arguments in the previous section, one can have 

l^n- nl^i <A(i,|/ l | Y ^ +a + ||V«|| Yja ^ +1 ) 

Then (15.91 follows from (15.191 . the Sobolev interpolation and Young’s inequality. 


□ 


6. Conormal estimates 

We shall derive a priori estimates on a time interval [0, T £,a ] on which it is assumed that 

( 6 . 1 ) 


dz<p > T 7 , \h\ 2,00 < — and g - dfq > ^ on {z = 0}. 
2 c 0 2 


Note in particular that this will allow one to use Lemmas 1C.II and 1C.21 

To derive the higher order energy estimates, we shall use the good unknown V a = Z a v — 
dfvZ a i|, a / 0. A key point is that the control of V a and Z a h will yield a control of Z a v: 


\Z a v\\ < 


Define 


\V a \\ + A ||Vc|| L oo N ) I Z a hU , \\V a \\ < \\Z a v\\ + A ||Vu|| Loo ') \Z a h\_i . 

V c o / 2 V c o / 2 

( 6 . 2 ) 


Aoo(i) = A ( —, |/l(t)| x -+5 + ||u(t)|| Y ^+5 + P^WIIx^+4 + \\d z v(t)\\ v rn +2 + £2 


Co 


xt 1 


'zzv\'>)\ | L°° 


(6.3) 

It will be shown then that those functions A(-, •) defined in the previous three sections can be 
bounded by Aqo for sufficiently large m, and also the elliptic estimates of q and the smoothing 
estimates of h will be restated along the way. First, taking k = ^ for m > 6 in the estimate 
(I4.15P yields 

Ikllyf + ll V, ?ll Y ^ + \\dzz<l\\ Y ^-l (6.4) 

^ A , \h\ Y ™+4 + \h\ x rn +i + |M| Y m + 5 + 11 V "D 11 y ^ 3 + |M| x ^+4 + 11 ^ V 11 X ^ +3 ^ < ^-OO , 

while taking k = m — 1 in the estimate (|4.14p gives 

IMIx ™- -1 HVgllx ^- 1 T ll^^^llx "* -2 

— A ( —, \ h\^ r m+ 5 + \h\^rn +3 + ]|u|| v ^ + l + ||Vc||^m + 3 + IMI^+3 + 11 11 ^ +2 


(6.5) 


c 0 Y- 


1 Y"2" 


XT 1 


x (\h\ xm _i + cr |/i| X m— 1,2 + |M| X m + ||Vc|| Xm _i + e |^| xl „_ 1; 3 + £\h\^ m _ h 3 ^ 

m.-\ + C Nx™- 1 - 2 + IMIx™ + ll^llx™- 1 + £ l U l x m-l,3 + £ |^| xm _l,3 ^ . 


< A 0 

Here one has required m > 6 so that by Sobolev’s inequality 

\h\ m _i_4 < \h\ ra+5 < Ihl m 
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and that by the anisotropic Sobolev embedding estimate (1A.5D 

|MI y -+ 3 +1 £ IMI X l+5 + \\d z v\\ x in +4 . 

It can be checked easily that all the functions A(-,-) defined in Propositions 15.ll and HOI and 
Lemmas 13.3113.51 and 13.61 are bounded by Aoo, due to (16.41) . Moreover, Proposition 15.11 implies 
that 

e\h(t)\ 2 xrn _ 1: s <e\h(0)\ 2 xm _ lt s Aoo^H^j+el/il^^) (6.6) 


and 


<£|M0)|^,1 +J o Acx> (e 


U|(L i + £ |/i| 2 . m 1 


. , _.... . , . (6.7) 

Proposition 15.21 together with (16.51) and the trace estimate |<?| xm _ lj i < ||g|| Xm _i + ||Vg|| Xm _i, 
yields that 


2 il|2 


O’ \ h \ 


<A n 


-1 + ^ l^"lx m_1 ’ 2 + IMIx™ + ll^^llx" 1 - 1 + £l Him—1,3 + £l N 


2 I „, 12 


.2 | u 1 2 


( 6 . 8 ) 


2 i ± _ 

Since \h\^ m -i ,2 < \h\ 5 m _ 1 5 , one may improve (16.81) by using Young’s inequality to get 

X ’ ^ 

C 2 |^-| 2 m _i,5 <Aoo (|/i| 2 m> _i + IMIxm + ll^z^llx " 1 - 1 + g2 + e 2 |^| 2 m — i,3 ) . (6.9) 

Note that in the following we will use frequently these L°° bounds involved in A^. 

6.1. Basic L 2 estimate. We start with the estimates of (v,h) itself, that is, the case a = 0. 
Proposition 6.1. For any smooth solution of (USD, it holds that 


\v{t)\\ 2 + g\h(t)\ 2 0 + a\h{t)\l + £ / ||Vu||‘ 


( 6 . 10 ) 


< 


Ao (\\vo \\ 2 + \ho\l + a\h 0 \l + J |M| 2 ^ 


Proof. Standard energy identity yields 
1 d 


2 dt 


[ \v\ 2 dV t + e [ \S^v\ 2 dV t = / (2eS lp v - ql) 
Jfl JQ J z =0 


— qd) N • v dy — 


f (2 eS‘ p v — ql) e% ■ v dy. 

J z=—b 

( 6 . 11 ) 


The Navier slip boundary condition implies that 

— j (2eS‘ p v — ql) ■ v dy = — j 2e{S v ve^)iVi dy = — 2 ke / \v\ 2 dy. 

J z=—b J z=—b J z=—b 

While the dynamic boundary condition and the kinematic boundary condition give 


( 6 . 12 ) 


/ {2eS v v — ql) N • vdy = — / (gh — aH)N ■ v dy = — / (gh-aH)dthdy (6.13) 
I z =0 J z =0 J z =0 

= -\j t f z _ o 9 \ h \ 2 + (V 1 + |Vh| 2 - l) dy. 


Consequently, 


7^ ^ + / 9 ^l 2 + 2(7 (V 1 + l V/l l 2 ” X ) dy ) ( 6 -!4) 

+ 2 e f {S^vl 2 dVt + 2ke f \v\ 2 dy = 0. 

J Q J z=—b 


VT+WW- 1 > 


1 1 




|v/i | 5 


Note that 
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due to (16.11) . and the trace estimate 

M Li({z=-b}) ll Vv ll INI + IMI 2 - 

Hence, (16.101) follows from ()6.141) . Lemma 1C.21 and Cauchy’s inequality. □ 

6.2. Estimate of ( Z a v , Z a h) for cto < m — 1. Next, we derive the energy estimates of 
( Z a v, Z a h) for 1 < \a\ < m and ao < m — 1, that is, except the cases a = 0 or «o = m. 

Proposition 6.2. Any smooth solution of (11.91) satisfies the estimate 

IM*)llxm-l.l + IM^Ix™” 1 . 1 + ° IMOIx™” 1 . 2 + £ + j Q £ + ° 2 l^l^m-1,5 

< A 0 ^||u(0)||| m _i,i + |/i(0)|| m _i,i + a |/i(0)|xm-i. 2 + £ |/i(0)|^ m , 


(6.15) 


/ 

Jo 


+ / A c 


i”-U + \h\tm-l + a Nx-A + e \h\^-l 3 + Ibllxm + 


"7 1 1 'X 

Proof. The energy identity for the equations (j3. 131) (13.141) yields 

~ [ \V a \ 2 dV t + 2e [ | S v V a[2 
2 at Jq Jfi 


V° | 2 dV t + 2e / \S*V a | z = X 0 Q + X fe Q + 7eg + ft' s , 


where 


T“ - 

x o — 


[ (2£S*V a - Q a l) N -V a dy 

Jz= 0 

[ (2 eS*V a - Q a I) e 3 • dy, 

J z=—b 

K%= [ (( d‘£vV‘ p vZ a ri-C a (T)-C a (q))-V a -C a (d)Q a ) dV t , 

Jn 

n a s = [ (eV a (S^v) + • ( £ a {v )) • V a dV t . 

Jn 

We first estimate Xfi. The boundary condition (13.391) implies 

[ (2 eS^V* - Z a ql + dfqZ a gl) N ■ V a dy 

Jz= 0 


z? = - 


/ ~rGL _ 

X 0 — 


(6.16) 

(6.17) 

(6.18) 

(6.19) 

( 6 . 20 ) 

( 6 . 21 ) 


[ —(g — dfq)Z a KN ■ V a dy + [ oZ a HN ■ V a dy + [ eC a (B £ )-V a dy 

J z =0 J z =0 J z =0 


L 


■V a dy— [ 2 eZ a h df (S v v) N • V a dy. 

Jz =0 


2eS‘ fi vUZ a N 

/ z=0 Jz =0 

By (13.401) . the third term in the right hand side of (16.211) can be bounded by 


eC a (B £ ) ■ V a dy 


I z =o 


< AqqS (|u|xm-l,l + |/l|xm-l,l) |H Q |o 


Due to (1A.11D . it holds that 


/ 

Jz =0 


2eS' ¥ TnZ Q N ■ V a dy 


Note that involves v^ll^zz^ll L°°, one has 


2eZ a hdf (S ip v)N-V a dy 


>z =0 


< 2e\Z a S7h\_i \S^vUV a \i 
2 2 

|^|.3 Hi- 


<2e\Z a h\ 0 \\d‘?(S*v)N\\ LOO \V a \ 0 


( 6 . 22 ) 

(6.23) 


(6.24) 


< A 00 £2|/ l | xm _ 1 , 1 \V a \ Q . 

For the first gravity term, one may use the boundary condition (13.471) to rewrite it as 

[ -(g-dfq)Z a hN-V a dy 
Jz=0 


(6.25) 
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I z =0 


-{g - dfq)Z a h ( d t Z a h + v y ■ V y Z a h + 8%v • N Z a h - C a {h )) dy. 


Integrating by parts in t and using (16.41) lead to 

1 d_ f 

2 dt J z= q 


[ -(g- dfq) z a hd t z a hdy <-~ [ (g - dfq) \Z a h\ 2 dy + . (6.26) 

J z =0 


A Aqq — L 1 


The integration by parts in y gives 

[ {g- dfq) Z a hv y ■ V y Z a h dy 

Jz= 0 

Due to (13.481) . one has 

[ (g- dfq) Z a h (dfv ■ N Z a h - C a (h )) dy 

Jz=0 

Hence, in light of the estimates (|6.26l) (16.281) . one may conclude from (16.251) that 

[ -(g-d*q)Z a hN-V a dy 
Jz= 0 

< -J id ~ dfq)\Z a h\ 2 dy + A 00 + Mx^-i) . 

To deal with the second term involving surface tension, one has by (15.131) that 

[ aZ a HN -V a dy = [ a\7 y 
Jz= 0 


(6.27) 


< Aoo \h\ (\h\ + M X m-1 ). (6.28) 


(6.29) 


' z=0 


. ( _ V y Z ° h _ Vyh • VyZ a h c(B a ) \ N -V Q dv 

{VTTWJF vATTW 3 ” ' J 

(6.30) 


where C(B°) is defined by (15.141) . By Lemma I3l6l one may deduce 

f <?Vy • C{Bf) N • V a dy < a |V y • C(Bf) | 0 \d t Z a h + v y ■ \7 y Z a h + d%v • N Z a h - C a (h)\ 0 

Jz =0 

< Aoo (7 \h\ Xm-1,2 (N X m,l + N Xm-1,2 + \h\ X m—1,1 + M X m-1 ) (6.31) 

< Aoofj \h\ X m —1,2 (H X m,l + M X m-1 ) • 

To study the other two terms, one rewrite it as, by using the boundary condition (13.471) again, 
( VyZ a h V yk ■ \/ y Z a h , 


[ 

Jz=0 


V y h N -V a dy 


t 3 v v 


<rV 


>z =0 


I _V_ _ jj _ » 

Wl + lVy^l 2 ^1+WfW 

( S7 y Z a h _ V y h • V y Z a h ^ 

' (w + |v„ft|2 ymwT 3 " J ( ‘ + "' “ ] y + Bi 


(6.32) 


where 


n * = 


cV,, 


( V y Z a h 


+ 


\/ y h ■ \7yZ a tl 


J 2=0 V V^+WfW Y / 1 + I Vy/ip 3 V / 

It follows from an integration by parts, (13.491) and (13.501) that 


V y h (dfv • N Z a h - C a {h )) dy. (6.33) 




X7 y Z a h 


V y h ■ V v Z a h r , 
—- 


V'i + IVy^l 2 ^l + |V/l| 2 


V y d>-NZ“h-C c 


(6.34) 


+ cr 


V„ 


( \7 v Z a h Vyh ■ V y Z a h 


I _V_ _ .v _V_ 

V\A + |Vy/l| 2 ^/1 + IV/ll 2 




\Z a ~ ai Vy ■ VyZ ai h\^ 


< Aoo (°" |h| X m-i,2 (|/i|xm-i,2 + M X m- 2 . 1 ) + fr|h| xm _i i 5 • 
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Integrating by parts in both y and t, one finds that 

f f X7 y Z a h Vyh ■ X7 y Z a h 

U w ' W! + lV y ^" > /rT|W 3VW 

( \S7 y Z a h \ 2 \V y h- V y Z a h\ 2 




where 




2 .L-, 


a a 


2=0 



y 7 ! + 


^i + lVyhl 2 ^/l + lV^I 2 


V„/i ] d t Z a hdy 

1“ ) dy + Kg 2 . 


(6.35) 


Vi + iv^l 5 


Vj,Z“/ir (6.36) 


Vyh-\7yZ ^ ydth . y yZ * h ] 


\/i + Iv.fiF 

— AooCT |/l| X m_l,2 ■ 

Similarly, the integration by parts twice yields 


crV, 


'2 = 0 


• ( /Tffvy \|2 “ ^7====J V ^l % ' dy < A oo<7 • (6.37) 

\ V 1 + |Vj/h| 2 ^/1 + lVhl 2 J 


Hence, by the estimates (16.3111 . (16.3411 (16.37(1 . one may conclude from (I6.30H that 


Id f ( \VyZ a h\ 2 \Vyh-VyZ a hf 


/ aZ°H~N ■ V a dy < -- / a . ___ _, 

/ 2=0 2 dt J z=0 ^ y/l+WyW y/TTWyW" 


dy (6.38) 


+ A 

OO (o-|h| ^rn —1,2 1^1^771,1 + C7 

Note also that Lemma lA.41 implies that 

< 


i-i.i M'S’-rn-l, A 


l^lo ~ |v|xm-l.l + Aoo and \ V °\l^ M x m-1,§ + A °° N x m —1,3 • 


(6.39) 


Consequently, plugging the estimates (16.2211 (16.2411 . (16.2911 and (16.3811 into (16.2111 . by (16.3911 and 
Cauchy’s inequality, one may finish the estimates of Iq as: 

( \VyZ a h\ 2 \Vyh-VyZ a h\ 2 ' 


T a < — 
x 0 — 




dy (6.40) 


yi + |V^| 2 y/l + \V y h\^ 

+ Aoo + |v| X m-l + O |/l| X m-l,2 N X m,l + <7 I ^ I x m—3., § M x m-1,£ 

+£ (y|x—W1 + l^l x m-l,3 + W x m-1,3 M x m-1,§ ) ) ■ 

It also follows from (13.421) and (|6.39jl that Z“ admits the following bound: 

2? = - [ (2eS^V Q e 3 ) i ■V l a dy = - [ 2e {kV« - S a (v) i3 ) ■ V? dy (6.41) 

J z=—b J z=—b 

< Aoo£ (|H“|o + |w| Xm _i,i + |h| xm _i,i) |I4 a | 0 

< Aoo£ (\v\ Xm-1,1 + W Xm-1,1 ^ • 

Next, the commutator Kc is estimated by using (13.18(1 . (|3.15l) . (13.161) . (16.2(1 . (|6.5p and (|6.4H 


as 


K a c < Aoo ((\\Z a r,\\ + ir (T)|| + \\C a (q)\\) ||H“|| + \\C a (d)\\ ||Q“||) (6.42) 

< Aoo (jh| xm _i + ||u|| X m-i + ||Vu|| Xm _i + ||Vg|| Xm -i) (iMIx™- 1 . 1 + 

+ Aoo (||Vv|| xm -i + |h| xm _i i j ^||g|| X m-i,i + |h| xm _i i ^ 

< Aoo ^|h| xm _l + a \h\ X m-l,2 + ||u || X m + ||Vu|| Xm -l + £ y| xm —! 3 + £ |/i| xm _ 1; 3 ^ 
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x (jMIx™" 1 ’ 1 + IIVpHx ™- 1 + ■ 

It remains to estimate the commutator TZg. First, it follows from the integration by parts, 
(13.1711 and (13.421) that 

[ eV^(£ a (v)) ■ V a dV t = - [ e£ a (v)- VV a dV t + [ e£ a {v)N • V a dy (6.43) 
JVt J£l J z =0 

< Aooe ((||V U || X _ 1 + ||VF“|| + |F“| 0 ) . 

Next, for the hrst term, one actually has to estimate 

n a Si = e [ C‘j{S <p v)i j Vj‘dVt (6.44) 

J n 

= £ [ C?i(S‘ p v)i j Vj*dV t + e [ Cf^SPv^jVfdVt + e [ C« 3 (S*v) l3 V«dV t 
J Q, J £1 J O 

:=K a s ’l + K a s ? + K a s f 

due to (|3.2I) . For 77.(h 1 , by (|3.3I) . it suffices to estimate terms like 

e [ Z^(^)(Z^d z ^v) ij )V J a dV t , 

Jn C'zf 

where (3 and 7 are such that (3 / 0, 7 / 0 and \j3\ + |y| = m. By using (13.lip , one can reduce 
the problem to the estimate of 

e [ c 1 ZP(^-)d z (Z^iS^v)^) V a dV t 

Jn VzV 

with (3 as before (thus |/3| < m — 1) and |y| < I 7 ) < m — 1. The integration by parts shows that 
it suffices to estimate three types of terms: 

h=e f ZP$V) Z' y (S <p v)i j d z V?dV t , 

Jn VzP 

X 2 = ej^ (dzZP^^Z^S^ijVfdVt, 

and 

l 3 = s [ ZP(^) Z 1 (S tp v)ijVj >l dy. 

J z =0 OzP 

For Z\ and X 2 , since f3 7 ^ 0, it follows from (IA.1I) . (13.71) and Lemma IB.21 that 

\h\ < Aooe (llVt;|| Xm _i + \h\ xm _ lt i) ||VF“|| 

and 

|2 2 |<Aoo£(||Vu|| Xm _ 1 + |h| xm _ 1 ,3)||y a ||. 

By (IA.1I) , /3 / 0 and Lemma 13.41 it holds that 


|231 ^ -Aqq£ (|h| xm _i,i + | V v | xm — 1 ) | F |q ^ Aqq£ (|h| X m— 1,1 T |u| X m— 1 , 1 ) | F |q 


Consequently, one can get from the previous three estimates that 

Vp|| X m-i + ||VF“|| + (||Vt ;|| Xra - 1 + ||F° 

+(i^-ix ™- 1 ’ 1 + i^ix m_1,1 ) in 0 ) ■ 


V al 
/C Si 


A AqqS 


The estimate of TZg’f is straightforward, one gets from the definition (|3.4I) that 


77' 


a, 2 
Si 


< Aooe" H ¥m -a \\V° 


To estimate 77)h 3 , one derives from (13.51) and (13.111) that 




< Aqo^Z |/l| _j 1 \\v c 


(6.45) 


(6.46) 



























28 


YANJIN WANG AND ZHOUPING XIN 


Note that one has used again in the previous two estimates the fact that A^ involves e 2 |<9 zz u|| ■ 

For the term 

e [ ^V a [Z a ,d z }(S^v)dVt, 

Jn dzV 

performing an integration by parts and using a similar arguments as for 1Z 1 Si show that 


dup 

d z y 


V a [Z a ,d z \(S^v)dV t 


<Aooe||Vt;|| xm -i(||V a || + ||Vy c 


Consequently, 

K a s f | < A^ (e-2 \h\ xm _^ ||F“|| + e HVuH^ (\\V a \\ + || W“ 
It then follows from (16.45f) (16.4711 that 

((^ £ l^l x ^-l,3 + £ || Vu|| xm - 


V a (S <fi v) ■ V a dV t 


< Ao 




(6.47) 


(6.48) 


+e (||Vu|| xm _i + |7i| xm _ 1 .1 ^ ||VI7 Q || + (\h\ xm -i,i + |u| xm _i,i) |V r “| 0 ^ . 
This, together with (16.431) . (16.21) and (|6.39l) . implies that 

T^s < A oo ((e* |/i| xm _i,i +s\h\ xm _ 1 '3 +e||Vu|| xm _ 1 ) (|M| xm _i,i + |h| xm _ 1: i) (6.49) 

+e (||Vu|| Xm _i + |7i| xrri _ 1> ^^ ||VV r “|| + |h| Xm _i,i + M X m-i,i) • 

We can now finish the proof of the proposition. By the estimates (I6.40D - (I6.42D and (16.491) . 

i_ 1 

the trace estimates < |M| X m-i + ||Vu|| xm _i and M xm -i,i < ||Vu||| m _ M ||u||| m _ 1:1 + 

||u|| xm _i,i, using Cauchy’s inequality, one may deduce from (|6. 161) that 


1 d_ 

2 dt 


£ a + 2e 


T Aqo 


[ \S^V a \ 2 dVt (6.50) 

J n 

X m-i,i + u |/i| xm _i, 2 |h| xm ,i + a\h\ _j 5 (||u|| xm _i,i + || Vu|| xm -i) 


+ + l|Vu||| m _i + ||u|| Xm + e |h|( 


+ 

+£ 


(|MIx"*-M + \\Vv\\ X m-l + |h| xm _^3) 
.-i + |h| xm _ 1 ,i)||VI/°||), 


where 
£ a : = 


V a \ 2 dV t + I (g — dfq) \Z a h\ z + cr 


a 7 12 


| VyZVtf* \Vyh-X7yZ a h\ 2 


6=0 


^V^+WyW ^TT\%W 


dy. (6.51) 


It follows from (16.21) . the Taylor sign condition in (16.11) and (15.181) that 


| Z a v\\ 2 + \Z a h\ 2 0 + cr \Z a h\l < A 


£ a . 


One can use the Korn inequality of Lemma 1C.21 and (16.21) to get that 


|W C 


< At - 

To 


dv t + |H 


\K 


On the other hand, by the definition of V a , 

£\\VZ a v\\ 2 < e||VKT| 2 + A c 


£Mhll 2 m _ 1 ,3 +" L " 2 


m — 1,4 


rm— l . rj - rsj 


Then integrating (16.501) in time, using the trace estimate |u| 
and Cauchy’s inequality, together with (16.101) . one deduces that 

||u(f)|| Xm _ 1 , 1 + |h(i)| xm _i,i + a I^WIx™- 1 . 2 + e [ l|Vu|| Xm _ M 


< 


l|Vu|b 


+ M 


(6.52) 
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< Ao ^||f(0)||xm-i,i + |^(0)lx m -i.i A a IM0)lx m -u 2 

+ L A °° + (T |/l| xm -l,2 + a l^l x m-l,| (ll' U llx ro - 1 > 1 A IIVpIIx-O 

+ 1 + ||Vu|||m-l + |M||m +£|^|^. m _ li 3^ • 

This, together with (16.61) . (16.91) and Cauchy’s inequality, leads to (16.151) . 


□ 


6.3. Estimate of (d]f l v,d'jf l h). We now derive the energy estimates of (d™v, d™h), that is, the 
case ao = m. 


Proposition 6.3. Any smooth solution of (11.91) satisfies the estimate 


3MS 


I (llM 2 + IWo + ^IW? + ^l^i ) 2 + / ( £ Jo im 

< tA o (n^(o)ii 2 + whml +* whml +^ iMo)i 2 m ,i ) 2 

+ t ( [ Aoo + CT |/l|xm,l + ]|w]|x"» + ll^z^llx” 1 - 1 + £ ll^llx™- 1 - 1 + £ 

2 m _l,l + ° |^||m-l,2 + ||u|||m-l,l + ||0 2 u||xm-2^ 

.-1 + a \h\x™,l + IMIx™ + ll^z^llx" 1 - 1 + £i + £2 1,§) ' 


(6.53) 


L 


+ / A 0 


Proof. In the current case, (|6.16l) can be restated as: 


1 d 


o , , \V m \ z dV t + 2s | S^V m \ z dV t =l™ +l™ + K 7 S + liy} + Kff, 

2 dt Jn Jq 


where 


q-m _ 
x 0 — 


[ (2 eS v V m - Q m l) N • V m dy , 

Jz= 0 

[ (2 eS^V m - Q m I) e 3 • V m dy, 

Jz =0 

[ C m {dWqdV t , 

Jn 

TVS = [ ((dfv • V^vd^ri - C m (T) - C m {q )) • V m + C m (d)dfqdf n r ]) dV t , 

Jn 

TVg= [ (£V m {S tp v)+£W‘ p ■£ m (v))-V m dV t . 


q-m _ 

H ~ ~ 


K% = ~ 


(6.54) 

(6.55) 

(6.56) 

(6.57) 

(6.58) 

(6.59) 


Here V m = d]fiv — dfvdffi'q, Q m = d™q — dfvd^r], and C m {-) are those commutators C a (-) for 
the case ao = m. Note that we have singled out the term 7Zq from TZff. 

We first estimate Xff , which can be rewritten as (similar to (16.211) 1 


q-m _ 
x 0 — 


[ -(< g-dfq)d] n KN-V m dy+ [ adf 1 HN ■ V m dy + [ eC m (B e ) ■ V m dy (6.60) 
J z =0 J z =0 J z =0 

- [ 2eS‘ fi vIld™N ■ V m dy — [ 2 ed™h d'f (S^v) N • V m dy. 

J z=0 J z =0 


Following the analysis in (16.221) (16. 24p . one can bound the last three terms in (16.601) by 


Aoo (e(Mx™-i.i + Wx™-m) \V m \o + e l^l x m A l^ m l| + £2 l^lx m \V m \o) ■ 

As (16.291) . one deduces 

[ -(g-dfq)drhN-V m dy 

Jz =0 


(6.61) 


(6.62) 
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1 d 


<—2 J t J z=o ($ - d tq) \ d T h \ dy + A c 


2 , | .2 

I ^hrm-1 


However, as explained in Section [21 one can not use the arguments leading to (16.381) to estimate 
f _ 0 ad™HN ■ V m dy since there is one half regularity loss for d™h so that it is difficult to 
control the following term, after using the kinematic boundary condition, 


ad™Hmd t N-8 t 


\7TL — 1 


v dy. 


(6.63) 


' 2=0 


The crucial observation here is that, this term will be cancelled out from estimating the term 
1Zq defined by (16.571) . So the estimates of this term will be postponed till we estimate 7 Zq. 
Similarly as (16.41|) . IV 1 admits the bound 

T™ < AqqE (|H m | 0 + |u| Xm _i,i + |h| xm _i,i) |H m | 0 . (6.64) 

We now estimate TVq- Note carefully that there is no any estimates of d™q, so one needs to 
integrate by parts in t. To continue, one needs more explicit expression of C m (d). Indeed, we 
will use a variant of (13.21) . It follows from the divergence free condition that 

d z tp (d\vi + <92^2) + d z v • N = 0. 

Applying <9™ to the above and using the definition of C m (d), one gets that 

d z cpC m (d) = [<9™, N, -8 z v\ + [d™,d z r), 8 x v x + <9 2 U2] • (6.65) 

Moreover, to integrate by parts in t, one needs to single out in C m (d) the highest m — 1 order 
time derivatives terms and use the following decomposition 


with 


d z v>c m (d) = c m (d)x + c m (d) 2 + c m (d) 3 + c m (du + c m (d) 5 

C m {d)x =m<9 f N ■ d?~ l d z v, C m (d ) 2 = rndtd.rid™-'(d x v x + d 2 v 2 ), 
C m (d ) 3 =m3™ _1 N • 8 t 8 z v, C m (d )4 = d z r/d t (d x v x + <9 2 U2), 

m —2 

C m (d) 5 = £<& (afN • dr £ d z v + dfd 2 ri • + 

1 =2 


Accordingly, 


k % = - / (r(rf)! + r(d) 2 + r(d) 3 +r(d) 4 +r(d) 5 )ar9^. 

Jn 

The Hfth term in (16.701) can be easily treated by the integration by parts in t as 


- [ C m {d) 3 8Tqdydz = -j f f C m {d) b 8™- l qdydz + KZ 


with 


1 - 1,4 + 


Urm—2 


I ar\\ 


K?= [ d t C m {d) 5 d? l - 1 qdydz< A 0 

Jn 

Integrate by parts in t to write the fourth term as 

- [ C m (d) i 8] n qdydz = -^ r [ C m {d) 4 d] n - 1 qdydz + n £\ 

Jn dt Jn 

where, by further integrating by parts in z and the trace theory, 

TZ r 4 n = f m (8™d z r]dt(divi + d 2 v 2 ) + 8™~ l 8 z rj<% {divi + d 2 v 2 )) 8 r t n ~ 1 q dydz 
Jn 


( 6 . 66 ) 

(6.67) 

( 6 . 68 ) 

(6.69) 


(6.70) 

(6.71) 

(6.72) 

(6.73) 


(6.74) 


[ md™hdt {dxvx + d 2 v 2 ) 8™~ l q dy - [ md^rjd, {{8 lV 1 + 8 2 v 2 ) cf^) dydz 
J z=0 J 

+ [ md^~ l d z yd^{8ivi + d 2 v 2 )d™~ 1 qdydz 

Jn 
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< a- (iar/»i_i |ar _1 ?|i + ii« i|ar _1 ?Li + ||<9r 1 <MI i^r^n) 


< Aoo|/i| xm ,_i \\d? g|| Hl . 

Similarly, integrate by parts in both t and y to bound the second and third terms by 

- [ (C m (d) 2 +C m (d) 3 )d?gdydz = -j f (C m (d) 2 + C m (d) 3 ) d?~ x q dydz + 77^, (6.75) 


where 


V m - 
/c 2,3 — 


[ dt (C m (d) 2 + C m (d ) 3 ) 

w -!,i + IMIx- + livwllxa) liar 1 ?! 


(6.76) 


^ Ago 




Finally, we turn to the most delicate term, the one involving C m (d)i in (16.7011 . Integrate by 
parts in z first to get 

- [ C m (d)id™qdydz = - [ md t N • d?~ l vd?qdy + [ md z (d?qd t N) ■ d?~ r v dydz. (6.77) 
Jfl J z =0 J £1 

Then integrate by parts in t to obtain 

[ rnd z (d?qd t N) • d?^v dydz = ^ [ md z (df V t N) • d?~ l v dydz + 11?, (6.78) 

Jn dt Jn 

with 


K[ n = - f rnd z ffi-'qdt N) • d?v + md z {d?~ l qd^ N) • d?~ l v dydz (6.79) 

Jci 

< Aoo 11 n 11 x m IIW^II • 

Note carefully that we integrate by parts in z first rather than in t since there is no estimates 
of d?v on the boundary. This also indicates the difficulty in controlling the first term in the 
right hand side of (16.771) since one can no longer integrate by parts in t. Recall here that there 
was also one term out of control, that is, (16.6311 . Our crucial observation is that there is a 
cancelation between them since q = gh — aH + leS^vn • n on {z = 0}. This motivates us to 
estimate together the first term in (|6.77l) and the second surface tension term in (16.601) . by the 
kinematic boundary condition, 


[ ad?HN ■ V m dy - [ md t N ■ d?~ 1 vd?qdy 
J z= 0 J z= 0 

= / ad?H (N • V m + md t N • d^v) dy 
Jz=0 


/ . 
Jz =0 


m&N • d?~ l v ( gd?h + ‘led? (S v v n • n)) dy 


ad?H 


I z =0 


+ Vy ■ X7yd?h + dfv • N d?h - C 


/ ■ 

Jz =0 


mdtN ■ d? l v ( gd?h + 2 ed? (S v v n • n)) dy , 


where C m (h) is the commutator C a (h ) defined by (13.5311 for the case «o 
last term in (16.8011 can be estimated as follows, thanks to Lemma 13.41 


(6.80) 


(h)) dy 

= m. Note that the 


m<9<N • d? l v ( gd?h + 2 ed? (S v v n • n)) dy 


lz =0 


< A 00 |a£ n - 1 u|i (\d?h\_i + e\d?(S*vn- n)|_i) 

< Aoopr-'H, (lafftl-j + £| Vi - 


(6.81) 
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The integration by parts and (13.501) yield 


' 2=0 


< a 


ad™H dfv ■ N d™h - C m 


dy 


(6.82) 




{- 


S7yh 


Vy dfv ■ Nd^h - C m 


< AooCJ |/l| X m,l (\d™h\i + |/l|^m-l,2 + |f| X m-2,l) ■ 

It follows from (15.131) and (I5.14|) that 

J aZ a H ( d™ +1 h + v y • V y d?h) dy (6.83) 

= L ^ ■ (t™? - + C(K) ) (8r+lfc+ *■ ^ h) *■ 

where 

1 


C(B™) = - 


dl 


,m —1 
t > 


Vyh 


71 T\%W 

Similarly as (16.351) (16.371) . one can deduce that 


v y d t hv v h+ dr, 


Vi+WW 


,Vyh 




( M2 - (AT 1 * + % . v/r*) iy 


2=0 \ V^M-WyW \J 1 + |Vy/l| 2 

1 d f z'^ + ^ 


(6.84) 


(6.85) 


< -- 


2 dt J z=0 y y/l + \V y h\ 2 ^l + \\7yh\ 2i 

Integrate by parts in both t and y to have 

[ aV y ■ (C(B?)) dT +l hdy < -4 f °C(Br) ■ Vyd^hdy + [ ad t C(B\ ?) • Vyd^hdy 

Jz =0 at Jz =0 Jz =0 

< -4 [ ^(B'a ) • V y drhdy + Aootr • (6.86) 

42=0 

One easily has 

[ a\7 y ■ (C(H™)) v y ■ Vydrhdy < A^a\h\ Xm -i, 2 \h\ X m,i ■ (6.87) 


' 2 = 0 


Hence, by the estimates (16.811) . (16.821) . (|6.85[) — (|6.87|) . one may conclude from (16.801) that 


J adrHN-V m dy- [ md t N • d^vd^q dy 


( 6 . 88 ) 


/ 2 = 0 


' 2=0 


< -- 


1 d 


( _ IV.'-W) dy _ * [ aCm . Vy8rhiy 


2 dt J z=0 y yj\ + | V y h\ 2 ^/l + | V y h\ 2,i J dtj z=0 

+ Aqo ^ct \h\^m,i + a \h\ X m,i |u| X m- 2 ,i + | d^ 1 v\i (\dl n h\_i + e |/i| xm> i + e 




This in particular finishes the estimates of the second surface tension term in (16.601) and 7 Zq, 
which can be stated as follows: 


[ ad^HN ■V m dy + K% 

Jz =0 

( \V y drh\ 2 \Vyh-Vydrh\ 2 ' 


(6.89) 


ld 


d *-Tt gm 


2dtj z , 0 \yrnwT yi + iVj/ip 

+ A-oo (a \h\ X m,i + a\h\ Xm ,i |u| Xm _ 2 ,i + 1<9™ 1 v\i + e |/i| xm ,i + £ 
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where 


g m = 


[ ( C m (d) 2 + C m (d) 3 + C m (d) 4 + C m (d ) 5 ) d™~\ dydz 

Jn 

- [ rnd z (d™- l qd t N) ■ d^vdydz + [ ■ Vyd^hdy. 

Jn Jz =o 


(6.90) 


It remains to estimates the commutators TZ™ and 1Z™. It follows from (I3.18|) . (13.1511 . (13.161) . 
(ESI) and (|6.4[) that 


TVS < Aoo ((IIM + \\C m (T)\\ + \\C m m • + \\C m (d)\\ IIM) ( 6 - 91 ) 

< Aoo + IMIx™- 1 + IIV^II^- 1 + ||V(?|| xm _i) (||F"'|| + ll^ll) 

< Aoo (N x m,-1 + G \h\ X m-l,2 + IMIx"-* + llV^Hxm-l + £ |t'| xm _l i 3 + £ |/l| xm _j,^ 

X (|MIX- + • 

Similarly as (16.491) . it holds that 

IZ 7 g <Aoo ^£ 2 |/i| xm _i,i + e |/i| xm _i,3 + £ ||Vu|| Xm _i^ ^||v|| X m + |/i| xm ,-i ) (6.92) 

+ £ ^||Vu|| xm _l + |/i| xrrl _ 1 1 ^ ||VT4 m || + (|/l| X m-l.l + |u| X m-l,l) (\h\ %m + Mx™)) ■ 

We can now finish the proof of the proposition. As a consequence of (16.611) . (16.621) . (|6.64l) . 
(I6.89p . (16.91ft . (16.921) and Cauchy’s inequality, one may deduce from (|6.16l) that 

-4 -£ m + — Q m + 2e [ \S^V m \ 2 dVt (6.93) 

2 dt dt J n 

< Aoo ^1 ^lx m + a Nx™a + IIVpII^- 1 + 11p 11x m + £ l^lx m + £ \h \ xm ,i 

+£ M x ^ (iMIx™ + HVullxa + 1^,1) + £ (l|Vu||xm-i + \h\ xm . 1>h ) HVW^Il) , 

where 


- [\V m \ 2 dV,+ / - V,) + a - M) *. (6.94) 

Jn Jz= o V V l + \ v v h \ y/1 + |V :y /i| 2 J 

Similarly as (|6.52|) , by the trace estimates 

i i 

M x ™,i % IMIx™ + HVullx- and |u| xm < ||Vu||£ m |M|f m + ||u|| xm , 
using Cauchy’s inequality and (16.71) . one can then deduce that 

ii«r»wii 2 + \<rm 2 o+° i owi?+= \m 2 r ,. +?f‘ nva? 


rjrTT, „. 112 


(6.95) 


< Ao (ll^(O)ll 2 + |«0)lo + * I d?h(0)\l + £ \h(0)\l 


~G r 


L 


+ / A c 


Note that 


~G m < A 0 


+ a Hx ™. 1 + IMIx™ + ll^z^llx™-! + £ ||Vu|| X m-l,l + £ l^l x m,^) • 


i-i ,\ + Mx™- 1 ' 1 + ll^ 1 ’llx m - 2 ) ||ar + G l^lx m -L 2 |^| X m,i) • (6.96) 


In contrast to the previous case, the difficulty here is that ||<9 t m 1 i/|| i/ i and hence — G m are not 
in L°°([0, T]) but only in A 2 ([0,T]). Our basic idea is to integrate in time twice. Indeed, we 
take the square and then integrate in time to have, by Cauchy’s inequality, 



(6.97) 
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< *a„ (iieromii 2 + \arm\l+ uarMo)i(+^ wo)i 2 „,.) 2 + J‘ is m i 2 
+ t (J hoe (|A||m +<r|A|xm,l + I|f||| m + ||O s w||xm-l + £ || Vf|| | m _ 1,1 + £|A| 2 m ,l)) 


It follows from (16.961) and (16.51) that 


r \Q m \ 2 < r 

Jo Jo 


Ar. 


+ Mx”*-M + II<%t||x ™- 2 


i,-l + a Nx ™- 1 . 2 (6.98) 


< 


+ IMIx™ + ll^^llx” 1 " 1 + £ l U l x m-l,| + e |/l| xrra _ 1 , 3 ^ l^lxm-1,2 |h||m,l^ 

j Aoo 1 + CT |/l|xm-l ,2 + Mxm-1,1 + ll^z^llx’"- 2 ) 

x (jh| 2 m _i + a + |M|| m + ||^u||xm-i + e 2 |t’| 2 m _ li 3 + e 2 | 3^ . 


We thus conclude (16.15[) by plugging the estimate (16.981) into (16.971) . 


□ 


7. Normal derivative estimates 

In view of the conormal estimates in Propositions 16.21 and 16.31 in Section [ 6 ] the next main 
step is to estimate ||(9 2 u|| Xm _i. 

Recall the definition of A^ from (16.31) and all the facts of the L°° controls elaborated in 
the beginning of Section [HI Note that Aoo involves only y/e\\d zz v\\ Lao . For the case without 
surface tension [331 that involves only the spatial derivatives, this is sufficient for deriving 
the normal derivative estimates since in such situation applying the product or commutator 
estimates to control the commutators resulting from the viscosity term needs only the control 
of \/~£\\d zz v\\ Lxl - However, in the current case that involves the time derivatives, following 
the arguments of [33] would require the control of y/e \\d zz v\\ Y k for some k > 1. Recall from 
Proposition 9.8 in [33] that deriving the bound of yfe ||<9 22 u|| i0 o requires a crucial use of the heat 
kernel and the first order compatibility condition <Sn|t=o = 0 on the boundary. Hence, to control 
yfe ||9 22 u|| Y fc, it seems to involve much more delicate use of the various properties of the heat 
kernel for the time differentiated problems; furthermore, it should require more compatibility 
conditions of initial data. 

Our key observation here is that since in the vicinity of the boundary the solution behaves as 
v(t,x) r^i v°(t,x) + y/eU(t,y,z/s /£), it indicates that there may be better control of ed zz v (and 
even ed zzz v\) in Sobolev conormal spaces. This is indeed the case as shown in the following 
lemma. 

Lemma 7.1. It holds that 

£\\ d zz v \\ Y ! n +1 + £ \\ d zzz v\\ Y rn <Aoo (7.1) 

and for k < m — 1 : 

£ ll < 9zz , y|| X fc < Aoo (e (ll^ullxM + N x fc,3) + ||<9w||xfc + IMIx*+i + W x h-i,-£ + ||Vg|| X fc) (7.2) 
and 

£ \\dzzzv\\^k-i < Aqo ^|| Vu||xfc-i ,2 + \h\^ k _ 1 s ^ + \h\^k,\ ( 7 - 3 ) 

+ ll^llx* + IMIxvi + l|Vg|| X fc-i,i + RVg|| xfc -i )• 

Proof. It follows from the first equation in (11.91) . (13.211) and (14.61) that 

ed zz v = I -e ^2 d z (E 3ij djv) - £ ^ di ( Eijdjv ) + d z ip (d t v + v y ■ A7 y v) + v z d z v + d z ipA7 v q 
33 V -?< 3 *<3,1 

(7.4)' 
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This implies that, since v z = 0 on the boundary, 


£ \\d zz v\\ Y ™ +1 < A (J-, £ ^||VV y n|| Y m + i + |^| y ^ +3 ) + \h\ 

l r+ 2 + \\^Q\\ 


M- 2 


(7.5) 


+ ||a 2 n|| Y ^+i + ||u| 


¥" 2 ^ 


YT 1 


< Aoo. 


+1 I A J1 -oo 


Applying d z to (17.41) . using the estimate (17.51) . one deduces 


e\\d zzz v\\ ™ < A — ,e 
To 


1 


+ V 2 V„u 


+ \hl 


iYij-+ 2 

+ \\d z v\\ Y ^+i + IMI Y ^+ 2 + ||<9 2 Vg|| Y - + ||Vg|| Y -^ < A c 

where one has used the fact that, since v z = 0 on the boundary, 


(7.6) 


z(y z d z v)\\ Y rn < \\d z V z d z v\\ Y rn + 11 U 2 <9 22 U 11 y ™ = \\d z V z d z v\\ Y rn + 
z y z || Y f- ||<9 2 n|| Y ^ + ||(9 2 n 2 || Y ^ \\d z v\\ Y ™ +1 . 


1 


z(z + b) 


v z Z 3 d z v 




(7.7) 


Combining (17.51) and (17.71) proves the estimates (17.11) . 

The estimates (17.21) (17.31) follow by a similar argument as that for (17.11) . Indeed, it follows 
from (17.41) and dAU) that for k > m — 1 , 

£ ||<9 22 u|| X fc < Aqo (jlVVyul^fc + \h\^ kt ^j + |/i| xfci i + ||3 z n|| X fc + |M| X fc+i + |/i| x *, +li _i + 

This yields (17.21) . Now applying Z % . i = 1, 2, 3 to (17.41) and then using (IA.1D again lead to 

£ ||$zz'P|lx fe - 1 . 1 — A-oo ^£ (lIVV^IIx*- 1 , 1 + |/l| xfc _ 1: S^ + |/fc| x fc_ 1> 3 (7.8) 

+ ll^z^llx*- 1 . 1 + IMIxM + \h\ x k,% + ll^lllx^-L 1 ) • 

Then applying d z to (17.41) and using the estimate (17.11) and the similar observation as in (17.71) 
give 


£ ||<9 222 u|| x fc-i < Aqo ^£ ^||V V J/ n|| xfc _ 1 + l^l x *-l,{j) "b l^l x fc-l, 

+ ||<9 2 n|| X fc + ||n|| X fe_i,i + \h\ xkt i + 

Combining (17.91) and (j 7.8 1) thus proves the estimate (17.31) . 

Lemma o then allows one to derive the normal derivative estimates. 

7.1. Estimate of ||< 9 2 u|| xm _ 2 . As in Section 8 of (33], one defines 

S n = II (S^vn — kx v ) . 

The main advantages of this quantity are that 

S n = 0 on {z = 0 , — 6 } 
and that the following estimates hold: 

Lemma 7.2. For every k > 0: 

||9 2 n|| x fc < Aqo ^||<S'n|| X fc + \h\ x k,l + IMIxM^ 

and 

||5 22 n|| X fe < Aqo ^||9 2 S n || X fe + |^| xfc| 3 + ||n|| X fc,2 


V<?|| X fc-i + ||Vg|| X fe_i ) . 


(7.9) 


□ 


(7.10) 

(7.11) 


(7.12) 

(7.13) 
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Proof. We start with the estimate (|7.12l) . The normal component of d z v is given by, due to the 
divergence free condition, 

r\ 

d z v-n = - n^j- (dwi + d 2 v 2 ) ■ (7.14) 

Then it follows from (1A.1D and Lemma IB. 21 that 

II 9 z v • n|| xfc < Aoo (||u|| x m + |/i| xfci i) • (7.15) 

It thus suffices to estimate the tangential components of d z v. Recall from the derivation of 
(13.371) that 

r\ 

Hd z v = — (pupXldiv + d 2 ipHd 2 v + 2S n — g^djV ■ nil d y , — «xllu) , (7-16) 

which follows from the derivation of (13.371) . Hence, by (IA.1D and Lemma IB. 21 one has 

lina^llxfc < Aoo (||v|| x m + \ h \ x k,l + Il'S’nllxfc + W d * v • n 11 ) • (7-17) 

Thus the estimate (17.121) follows from combining (17.151) and (17.171) . 

The estimate (17.131) follows from applying d z to ([7. 141) and (I7.16D . employing a similar argu¬ 
ment as for (|7.12l) and using (|7.12l) . □ 


In light of Lemma 17.21 we then turn to estimate S n instead of d z v. It follows from the first 


equation in (11.91) that 

dfV^v + (v ■ V*) + (V^u ) 2 + (V ^) 2 q - eA^v = 0. (7.18) 

Taking the symmetric part of (17.181) yields 

Sf STv + (v • V*) S*v + \ ((V^v) 2 + ((W)*) 2 ) + (V ^) 2 q - eAV(S*v) = 0. (7.19) 

Consequently, 

dfS n + (vV*)S n -eA*(S n ) = Fs (7.20) 

where 

F S = F' + Fl + Fl + F$ (7.21) 

with 

Fs = -\n ((V v u) 2 + ((V^u)*) 2 ) n + (d t U + V ■ V^n) S^v n + US^v (d t n + v • V^n) (7.22) 
+ /e(flf + v- V^)(xnu), 

F 2 = -n ((V ^) 2 g) n, (7.23) 

Fl = -£ (Am) s v v n - eUS^v A^n, (7.24) 

Fl = - 2edfUdf (S*v n) - 2eU (df (S^v) dfn) - £kA^{ x TIv). (7.25) 


Note that the pressure estimates in (16.51) indicate that at most, one can estimate only |,S n || X m -2 
and hence ||<9 z Hlx m - 2 this stage. However, we shall prove a control of ||9 z u|| Xm _i based on 
the vorticity equation in the next subsection. 


Proposition 7.3. Any smooth solution of (11.91) satisfies the estimate 


ISnH — 2 + £ / ||VS n || —2 

Jo 

rt 


(7.26) 


A Ao ||‘S' n (0)||xm-2 + J Aoo£- ^11 Vu|| X m_2,2 + \h\^ 

+ Aoo + CT 2 \h\ X m- 2.3 + 


I ^ 11 —IjI I — 2,1 + || An || Ym—2 



































\F s \\ xm - 2 < Aqo ||V(j , || Xm - 2 ,i, 

| A 1 ||| Xm _2 < AqoE | || Vv || xm _2 + |/ l | xm _ 2,5 ^ , 
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Proof. We start with the estimates of F$. It follows from (1A.1D and Lemma IB. 21 that 

K||*»- a — Aoo (||Vv|| X m-2 + |h| xm _i,l + ||u|| X m- 2 ^ , (7.27) 

(7.28) 

(7.29) 

l|-^s|| X m—2 < AooS [\h\^m-2 +2 ||V 2 f|| xm _ 2 + |A| xm _ 2 ,3 || V 2 'U'11 y m - 2 ^ (7.30) 

< Aoo (e ||Vn|| xm _2,i + ||5 2 r;|| xm _2 + |M| X m-i + + ||Vg|| xm _ 2 + W xm _ 2 ,3^ . 

Here in the second inequality of (I7.30|) . one has used (17.11) and (17.21) with k = m — 2. Hence, 

\\Fshm -2 < Aoo (e ( ||Vi;|| Xm - 2 ,l + H xm —2,5 ) + \\dzV llx"*- 2 + IMIx™- 1 + l^l x m-l,l + l|Vqi| X m- 2 ,l 

(7.31) 

(7.32) 

(7.33) 

(7.34) 


It follows from applying gC-Z a (d z ip-) for |a| < m — 2 to (17.201) . (I3.2ip and (14.61) that 
dfZ a S n + ('v ■ V?) Z a S n - eA Lp Z a S n = Z a F s + C s , 
where the commutator is given by 

Cs = Cg + C‘s 

with 


Cl = 


1 


^,8^ (dt + Vy ■ Vy) S n + [Z 0 , Vy] • VyS n + [Z“ , V Z \ 3 Z S n 


S d z ip 

:= Cl, + Cg 2 + Cg 3 


and 


cl = -£ 9 ^[Z“. V • (B y Sn) - £ v • ([Z“, E] vs n ) - V • (E[Z“, V] S„) 


(7.35) 


Cl, + C 2 2 + C| 3 . 

Since Z a S n = 0 on the boundary, one can obtain 

[ \z a s n \ 2 dVt+£ [ \V*Z a S n \ 2 dV t = [ ( Z a Fs + Cs)-Z a S n dV t . (7.36) 

z dt Jfi Jci Jq 

The right hand side of (|7.36l) can be estimated as follows. (I7.3ip implies immediately that 


/ 

Jn 


Z° 


(7.37) 


n to— 2 


*F S ■ Z a S n dV t < Aoo (e (ll^ullxm- 2,1 + ) + \\9 z v 

+ ||n|| xm _i + \h\^_ h i + ||V(j|| xm -2,i^ ||5, 

Next, we estimate the part involving Cg. (IA.2D yields 

II^Sill + < Aqo (|Mlx ™- 2 + II|lx m — 2 + l^l x m-2,i) • 

Additional care is needed to estimate ||Cg 3 ||, since d z S n can not be controlled. By expanding 
the commutator and using (13.111) . one sees clearly that it suffices to estimate terms of the form 

z^v z d z z-s r 

with |/3| + |y| < m — 2 , |y| < m — 3. Since 


(7.38) 


z^ Vz d z z^s n = 


1 


-Z?v z Z 3 Z^S n , 


z(z + b) 

which can be further rewritten as a sum of terms of the form 


1 


{z(z + b)' 


v z Z :i Z>S n , 


(7.39) 
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where are harmless bounded functions and \j3\ < |/3|. Indeed, this comes from the fact that 
Z 3 ^ ' If /5 = 0, since (13.221) implies that v z = 0 on the boundary, then 

1 


r~Z@ . _ 

^ V-(- + & ) 

If (3 / 0, one can use (1A.1D to obtain that 
1 


v z Z :i Z~>'S n 


\\d z v 


Z\\L° 


|S n | 


Z P Zl3 \z(z + b)' 

1 


v z Z 3 Z^S n 


< 


Z 


11 Z 3 Sn 


n-3 + 


^z(z + b) *) 

Observe that by again that v z = 0 on the boundary, 

1 


Z 


z(z + b ) 


||^3 | Sn||y m -3 • 


z 


z(z + b) 


m —3 ^ W d * V * 


On the other hand, since 

z ' 1 


,z(z + b) 
it suffices to estimate 

1 


< 

rs_/ 


1 


z(z + b) ZVz 


+ 


z(z + b ) 


Z^Zr 


+ b) 


Z p v 2 


z(z + b) 


< m — 3. 


Indeed, since v z = 0 on the boundary, it follows from the Hardy inequality that 


. 1 v .Z fi Zv z 

< 

d z Z f3 Zv z 

and 

, 1 . Z^v z 

< 

d z Z ,3 v z 

z(z + b) 




z{z + b) 




We have thus proven that 

|(f'S’311 — \\ d zV z \\ Y ^ ||5 n || X m-2 + ll-S'nll^ ll^'WzIlx™- 2 ■ 

By Lemma IB.21 it holds that 

II^Z^zIIy^T — A(X> Hearts || Xm _2 ^ Aqq ^ll^llx m_2 ffi ll^z'^llx™' -2 ~b • 

Hence, 


I ^-'S’3 11 — Aqo 


-2 + ||3.a|| xm -2 + \h\ _! 1 + II Sj 


n x m ~ 2 


(7.40) 


This together with (17.381) yields 

J C$ ■ Z a S n dVt < Aoo ^||u || Xm _ 2 + ||3 2 u|| Xm _2 + + Il'S'nllx" 1 - 2 ) ll^nllx ™- 2 • (7-41) 

Next, we shall estimate the term involving C|. As mentioned in the beginning of this section, 
we need to employ a different argument from Proposition 8.3 in [33]. Due to (13.111) to handle 
the term involving C| j, it suffices to estimate 

e [ Z^d z (EVS n ) ■ Z a S n dydz 

Jn 

with \/3\ <m — 3, which can be reduced to the estimate of 

e J (z p, d z EZ p "vS n + Z 0, EZP"d z VSn) • Z a S n dydz 

with (3 1 + f3" = /3. It follows from (1A.1I) . (17.11) and (17.21) with k = m — 3 that 


Z^'d z EZ^"VS r 


1 $ ll^zL'II m-3 £ ||V5 n || Xm _ 3 + || 9 2 Ll || Xm _3 £ || V5 n || m- 3 


(7.42) 


< Aoo ^£ ||5 2 u|| xm _3,l + ||9 2 u|| X m-3 + ||a|| X m-2 + l^m-2,-1 + ||Vg|| X m-3 + I ^ I —3, § ^ J 
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and (|7.3I) with k = m — 2 implies that 


zP'EZP"d z VS n \\ < \\E\\ y ^ e ||d 2 VS n || xm _ 3 + ||£|| xm -3 e ||^V5„|| yH ^3 (7.43) 

< Aoo (e ^||Vn || Xm _ 3 .2 + + |/i| xm _ 2 ,i 

-2 + IMIx m -2,i + ll^ ( ?llx m - 3 ’ 1 A ll^Vgllx — 3 ^ • 


+ WdzvW 


Thus, 


Cg\ • Z a S n dVt <Aoo (e (ll Vu|| X m - 3 ,2 + \h\ xm _ 3 ,s) + K 


(7.44) 


+ ll^z u ||x m - 2 E IMIx™- 2 , 1 + IIVgllx^-M + ||<9 2 V(/|| Xm _3 ^ • 


By expanding the second commutator C| 2 , one sees that it suffices to estimate terms of the 
form 


e / V • (V£Z 7 VS n ) • Z a S n dydz 


with /3 + 7 = a,/3/0. If 7 = 0 and hence /3 = a, since Z a S n = 0 on the boundary, one can 
then integrate by parts to get 


: f V • (Z a EX7S n ) ■ Z a S n dydz 
■hi 


s / Z a EVS n • VZ a S n dydz 
Jn 


(7.45) 


< £ \\Z a E\\ ||VS' n || ioo ||VZ Q S n || < Aoc \h\ 1 £2 ||VZ“5 n ||, 


where in the last inequality one has used the fact that Aqo involves yfe ||<9zzw|| i(X) . If 7 / 0 and 
hence 1 < |/3| < m — 3, then one can expand V-, by (IA.1D and Lemma mi to estimate 

e [ (z p \/EZ^VS n + Z p EZ^\/ 2 S n ) ■ Z a S n dydz. 


It follows from 1 A. If) , (17.111 and (17.21) with k = m — 3 that 

£ Z^VEZ^VSn < IIZV^II m -4 £llZV5„ll xm -4 + llZV£;il xm -4£|lZV5w.ll m -4 (7.46) 

< Aqo ^£ ||5 2 n|| Xm _ 3 ,i + ||5 2 n|| Xm _ 3 + |M| Xm _ 2 + |/i| xm _ 2 _i + ||V(7|| Xm _ 3 + |/i| xm _ 3j 3^ , 

and (USD with k = m — 2 that 


Z^EZ^V 2 S n 

< A„ 


+ d z v 


~ 11 Z E 11 ^ m ~ 4 £ ||ZV 2 S ' n || xm _ 4 + ||ZA , || Xm _4 £ ||ZV 2 iS l n | 

-3,2 + N x ™—3,§) + \h\ x m-2^ 

-2 + IMIx ™- 2 * 1 + IIV^IIx-M + ll^z^ f /llx m - 3 ) • 


(7.47) 


z u \\X rn ' 


Hence, 


/ 
J n 


Cs2 ' Z a S n dVt < Aqo (e (j|Vv|| X m-3,2 + |/i| xm _ 3 i s^ + |/i| xm _ 2 ,i + ||<9 2 


+ ||Vg|| Xm -3,i + ll^z^^llx ™- 3 j Il'S’nllx "*- 2 A ^-00 \h\ 
To handle the term involving C| 3 , due to (13.111) . one needs to estimate 


'O 11 ’J^rn — 2 + INI 2,1 
1 


1 £2 ||VZ Q S n | 
’2 11 111 


(7.48) 


: J V • ( EZ?d z S n ) • Z a S n dydz 


with \fd\ < m — 3, which can be bounded easily by using (17.21) with k = m — 3 and (17.31) with 
k = m — 2 so that 


/ 


C 2 S3 -Z a S n dVt < Aoo(e 


1-3,2 + N x m-3,§) + Wj 


■m- 2 ,i T Hoz^llx " 1 - 2 


(7.49) 
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+ IMIx^- 2. 1 + ||Vg|| Xm _ 3 ,i + ||(9 2 Vg|| Xm _3 ^ ||S' n || Xm _2 ■ 

In view of (I7.44j) . (17.480 . (I7.49p . one has actually proven that 

Jc 2 s ■ Z a S n dV t < Aoo (e ( ||Vu || xm _ 3 ,2 + |^l xm _ 3 ,§) 4- |^l xm _ 2 ,i + lldz^llx ™- 2 + IMIx ™- 2 . 1 

+ l|V g || Xm -3,l + ||^2|V(7|| Xm -3 ^ ||S n || Xm _2 + Aoo |^l x m-2,^ £2 l|V-Z’“iS' n | 


Consequently, by (17.371) . (I7.41|) and (17.501) . one deduces from (17.361) that 

~ [ \z a s n \ 2 dVt + £ [ \V*Z a S n \ 2 dVt 

* dt Jn Jft 


(7.50) 

(7.51) 


< A„o (£ (IIVu 


+ \hi 


J d - ||5 n ||3j- m _2 


+ \\d Z V Ym-2 + M 


+ |/i| x m-l,l + ||Vg|| Xm - 2 ,l + ||5 z Vg|| X m-3^ ||5 n || X m-2 + Aoo |^l x m- 2 ,^ £2 ||VZ“5 n ||. 

To conclude, one can use Lemma fC.ll to replace ||V ¥, 5 n || Xm _ 2 by ||V 5 n || Xm _2 in the left hand 
side and then use Cauchy’s inequality to absorb the last term. On the other hand, one can 
follow the derivation of (16.51) to get that 

||g|| X m- 2 ,i + ||Vg|| xm _ 2 ,i + ||<9 22 g|| xm _2 

< Aoo (j/l| xm _l,l + O l^lxm—2,3 + ||u|| X m-l,l + ||Vu|| Xm _ 2 ,l + £ || V U11 X m-2,2 + £ |^| xm _ 2 ,§) ) 

which reduces the order of time derivatives to the spatial derivatives. Finally, we integrate the 
resulting inequality in time to obtain (17.261) . □ 

7.2. Estimate of ||<9 2 u|| Xm _i. Note that Proposition 1 7.3 1 only provides the control of || 3 2 u|| Xm _ 2 , 
and this is due to the appearance of (V^) 2 ^ in the source term in the equation of S n . To get 
around this, a natural way is to use the vorticity instead of S n . 

Set to = x v. Then 


1 


u x n = - (V^v n — (V^u) 4 !!) = S v v n — (V^u^n, 


(7.52) 


and hence by (13.351) . 


lo x n = -8 n u — g l i(djV ■ n )d y i. 

Consequently, as in Lemma 17.21 one can get that 

H^uHxm-i < Aqo ^ll' y llx" 1-1 ’ 1 4" l^l x m-i,^ T IMIx" 1-1 ) • 

To estimate ||cu|| Xm _i, one notes that the first equation in (11.91) implies 

dfu + v- — u ■ = eAFui. 

As in the previous subsection, applying Q^Z a (d z tp-) for |a| < m — 1 to (17.541) yields 


Here F is given by 

where Cs is given, as in (17.33[) . by 

with 

1 


d?Z a uj + (v ■ V?) Z a oj - eAZ a u = F. 
F = Z a (yj ■ V^u) + C s , 


Cs -Cs+ Cs 


(7.53) 

(7.54) 

(7.55) 

(7.56) 

(7.57) 


4 = — [Z a ,8 z <p\ {d t + vyWy)u+ [Z a ,v y ] ■ VyU + [Z a , v z \ 8 Z UJ 
OzF 

:= C x sx + Cg 2 +Cg 3 


(7.58) 
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and 


C 2 S = -£^~[Z a , V] • (EVu) - £ 77 —V • ([Z Q , E\ Vw) - e-^-V • (E[Z a , V]w) (7.59) 

U z tp O z lf 0 z Lfi 

:= Cl, + C52 + C53. 

The main difficulty lies in that the vorticity does not vanish on the boundary. Thus, set 


where sloves 

' df<h + (« • ^Xh ~ ^Xh = F m 

< = 0 on {z = 0 , -b} 

Xh\t=o = w(0) 

and solves 

' + (u • - £A^w£ = 0 in Q 

< on {z = 0 , — 6 } 

XX = °- 

The estimate of u;"^ is given as follows. 


(7.60) 


(7.61) 


(7.62) 


Proposition 7.4. For \a\ <m — 1, the solution uj^ h of (17.6111 satisfies the estimate 


ll w n/i( i )l| 2 + £ / ||Vw“ ft || 2 <A 0 ||w(0)||| m —1 + [ A ac £ 2 


*- 2.2 + l ^| 2 m -2 , 5 ) (7.63) 


L 


+ A c 


+ OiS ¥m-l + U hfm-1,1 + 


.a 11 2 


Proof. Since = 0 on the boundary, it follows from (17.6111 that 

~ [ Kh\ 2 dV t + e f \VX h \ 2 dV t = [ F-^ h dV t . (7.64) 

We now estimate the right hand side of (|7.64|1 . (1A.1I1 implies that 

II Z (t O ■ V*V)|| < Aqo (||w|| X m_l + ||V^|| xm _l) < Aoo ^||w|| X m-l + ||Vv|| X m-l + l^l^m-1,^ ■ 

(7.65) 

To estimate the part involving Cs, one first can change S n into ui and m into m + 1 in (17.5011 to 
get 

j^Cl-^ dVt < Aoo (e (||Vu|| X m - 2 ,2 + \h\ xm _ 2 ,s) + lhl ^.,,1 (7.66) 

+ ||<9 2 u|| X m-l + ||v|| X m-l,l ) ||<V“J| + Aqo £2 ||V<J . 

For the part involving C(., one can change S n into u and m into m + 1 in (17.3811 to obtain 

Kill + Kill ^ (iMIx ™- 1 + IMIx ™- 1 + H x m— 1 , 1 ) • (7.67) 


The commutator C^., requires much more care. Indeed, one can not change m into m + 1 in 
(17.4011 since it would involve \h\ m 1 . As in (17.3911 . this commutator can be expanded into a sum 
of terms of the form 

csz " (^W«) z3ziu 

such that |/3| + |y| < m — 1, |y| < m — 2. As the arguments in previous sections, by (1A.1I1 and 
v z = 0 on the boundary that, one deduces 




z(z + b ) 


Z 3 Z' y uj 


<A r 


M 


+ 


Z 


z{z + b ) 


(7.68) 
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£ Aqo I W Ym-2,1 + 


1 


z(z + b) 

It follows from v z = 0 on the boundary and the Hardy inequality that, since v z = v ■ N — dtrj, 

1 


z(z + b) 


< 


^Z^z\\Xm—l 


(7.69) 


< Aoo (||^u|| xm -i + ||u|| X m-i + + ^2 \\v ■ d z Z a N - d z Z a d t T]\\ , 

\a\=m—l 

where the last term requires furthermore analysis. Due to (11.611 . it holds that 

\ z sv\ % \x(\Z\z)h\ 

where y has a slightly bigger support than y. Hence Z 3 acts as a zero order operator: 

l|V^ 3 r?||<Ni. (7.70) 

This yields that if a 3 ^ 0, one gains at least one derivative, and thus, 


• d z Z a 'N - d z Z a d t r,\\ < A^ 


a _2,3 + l<9^l x m-2.l) < A, 


OO • 


(7.71) 


Hence it suffices to estimate this term for the case with |a| = m — 1 and a 3 = 0. Note that 
(11.61) implies that 

1 ’ (7.72) 


(7.73) 


V= \ 1+ b) ^{z)* yh:= \ 1+ b) 

where * y stands for the convolution in the y variable and i/j is in L 1 (M 2 ). Then 
v ■ d z Z a N - d z Z a d t rj = - ^ (v y ■ ipz *y V y Z a h + ip z * y d t Z a h) 


- (l + 9 Z (Vy ■ lJj Z *y VyZ a h + l/j Z *y ~ T a = T a l + 7^2 • 


It is clear that 


I'TqlI || — A-oo 


,_!,1 + < Ac 


-h. ' 


(7.74) 


(7.75) 

(7.76) 


For the second term, one can separate it into two cases. For —b<z<~^, then 

Ta 2 = -J^Z 3 (v y ■ l/j Z * y V y Z a h + l/j Z *y S^/l) . 

It then follows from (17.701) that 

HT Q 2|| L 2({_ fe <,<_ |} ) < Aoo — Aoo ■ 

For — | < z < 0, T q2 can be written as 

Ta2 = - (1 + |) (v y (t, y, 0) • d z (^ *y VyZ a h) + d z (^ % d t Z a h )) + 77, 

where 

\K\<A 00 \z{z + b)\\d z {'<p z *Z a Vh)\<A 00 \Z 3 (ip z *Z a Vh)\. 

By using (17.701) again, one can get that 

ll^llL 2 (nn{-f <z<o}) A Aoo ■ 

To estimate the first term in (17.761) . one notes that 

Vy(t, y, 0) • d z (ip z *y V y Z a h) + d z {ij) z * y d t Z a h ) = d z (ip z * y (v y (t, y, 0) • \7 y Z a h + d t Z a h)) 

+ d z {v y {t, y, 0) - v y (t, y 1 , 0)) ■ ip z {y - y')\7 y Z a h(t , y')^ ■ 

For the second term in the right hand side of the above, one can employ the Taylor formula for 
v y to get that 


(7.77) 
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< Aoo 1/iL™-! 1 . (7.78) 


where i/j is still an L 1 function. This yields that 

sup dz([ (v y (t,y,0) - v y (t,y',0)) • fi z (y - y')V y Z a h(t,y') 

ze(~ |,o) ViR 2 

For the first term, one shall use (13.4711 to get 

v y ■ V y Z°h + d t Z a h = -dfv • N Z a h + C a (h) + V a • N, 

which implies, since |a| = m — 1, 

I Vy • VyZ a h + d t z a h\i < Aoo + \h\ x7n _^ 1 ) . 

It thus holds that 

\\d z (V’z *„ (y y (t,y, 0) • V y Z a h - d t z a h))\\ < I Vy ■ V y Z a h - 8 t Z a h | 1 (7.79) 

< Aoo + I^lx m- 1’0 ’ 


Collecting all the estimates (17.77H - (17.791) . one deduces from (17.7611 that 

II^II.L 2 (nn{z>-§}) - A °° (M x m-i,i + W x m-i,i) . (7.80) 

In view of the estimates (I7.68p . (17.6911 . (I7.71|) . (|7.74p . (I7.75P and ([7.8011 . by the trace estimates, 
one finally gets 


C 


S311 — A °o 


9 2 n|| xm _i + ||n| 


+ 


C _1 'i 


(7.81) 


As a consequence of (|7 .65 P . (I7.66P . (|7.67l) . (17.8111 and Lemma 1C.11 we deduce from (I7.64p that 
ll w h( i )l| 2 + £ J llVcv^H 2 < Ao ||w(0)|| Xm _i + / Aoo (£(||Vn|| xro _2,2 + |/ l | xm _2,5) (7.82) 

+ l^l x m-i.i + lld^llx™- 1 + IMIx™-M ) W^nhW + Aoo£2 |/l| xm _!,i IIV^II . 

We then conlude the proposition by using Cauchy’s inequality. □ 

It remains to estimate of (17.621) . Note that Lemma [3.41 and the trace estimate imply that 

y/s j |2I q cj|q < yfe / A. ^|n| xm _i,i + |/i| X m-i,i^ (7.83) 


< 


Jo 


Ar 


x— 1,1 L^ V’m—1,1 + M 


1,1 + |h|. 


Thus the main difficulty will be to handle the non-homogeneous boundary value problem, (I7.62p . 
whose boundary value is at a low level of regularity (it is Lf y and no more) due to (17.831) . This 
creates two difficulties: the first is that one cannot lift the boundary condition easily and 
perform a standard energy estimate; the second one is that due to the lower boundary estimate, 
one cannot expect an estimate of \\oj^\\ in L°°( 0 ,T) independent of e, as was well explained 
in Section 10.2 of [33]: in using the model of the heat equation, one may expect a control in 
H*(0,T) C L 4 (0,T). 

Proposition 7.5. Assume that A oc (t) + ^ < M for M > 0. Then for \a\ < m — 1, 

the solution of ([7.6211 satisfies the estimate 


hWl^L 2 — A(Af) / ^||n||| m _i,i + + e / ||Vn||j 


(7.84) 


Proof. Since the situation here is similar to Section 10.2 in [33], thus this proposition is a 
restatement of Proposition 10.4 in [33] . However, for completeness, we will sketch the main idea 
and some steps of the proof and state our estimates with a slight modification. 
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First, it is convenient to eliminate the convection term in (17.621) by considering Lagrangian 
coordinates. Define a parametrization of Q(t), X(t,-) : H i-> fl(t), by 

(d t X(t,x) = u{t,X(t, x)) = -) _1 oX) 

\*( 0 ,s) = $( 0 ,x), 

where the map <h(t, •) is defined by (11.41) . Denote J(t,x) = |det VX(i, s)| for the Jacobian of 
the change of coordinates. Then J(t,x ) = Jo(x) by the divergence free condition. Define 

Q a = o X), (7.86) 

where 7 > 0 is a large parameter to be chosen. Then solves 

ao(<9tD" + 7 !!“) — edi ( aijdjQ a ) = 0 in fl. (7.87) 

Here ao = | «/o I ^ and the matrix (a tJ ) is defined by 

(a t j) = | Jo 1 2 with Pij = diX ■ djX. 

Due to the assumption A 00 (t) + ||<9 22 u|| loo < M, Lemma 10.5 in [33] holds. Then the 


following estimates hold: 

II^11 w 1 ’ 00 11^ llw 1 ’ 00 — -^-O) (7.88) 

||VX|| LOO + \\d t VX\\ LOO < A (M)e tA ^ M \ (7.89) 

l|VX|| li00 + ||a t VX|| 1)00 < A {M)e^ M \ (7.90) 

IIV 2 X||i i0O + Vi||^ v2 X|| Loo < A(M)(! + t)e^ M \ (7.91) 

Indeed, the estimate (17.881) follows directly by J(t,x ) = Jq{x). Next, (17.851) implies that 

8 t VX = VuV$ _1 VX (7.92) 


and hence 

l|Vu|| L oo || VA'H^oo < A 0 + A(M) f ||VAT|| ioo . 

Jo 

This yields the first part of (17.891) by the Gronwall inequality. Next, applying one spatial 
conormal derivative to (17.921) . one can get 

l|Vu||i i00 ||VX|| liOD < A 0 + A(M) f ||VX || li00 . 

Jo 

and hence the first part of (17.901) follows from the Gronwall inequality. The estimates hold also 
for the time derivative in (17.891) and (17.901) by using again the equation (17.921) . To prove (17.911) , 
one applies \fed Xi to (17.921) to find that 

d t Ved Xi VX = sfe (3 Xi VuV$ _1 VX + V$ _1 VAT + VuW^VX) . (7.93) 

It follows from (17.901) that 

^||V 2 X|| loo < A(M)(l + f) + A(M) f ^\\^ 2 X\\ loo + A(M)e A ^ t f y/i ||V 2 u|| loo 

’ Jo ' Jo 

and hence, by using the assumption f* ^ < M, 

y/e ||V 2 X || < A(M)(1 + t)e A ( M ^ + A(M) f* y/i ||V 2 X || 1;00 

’ Jo 

and so the first part of (17.911) follows from the Gronwall inequality. For the second part of (17.911) , 
it follows by using again (17.931) and the fact that A^ involves the control of \fe || V 2 u|| ioo . 

As consequences of the estimates (I7.88|) (17.911) , one gets 

oo > m, 037 > m, (atj) > cold (7.94) 

by a suitable choice of m and Co depending on M for t E [0, 1] and that 

II (® 0 > ) || £,00 + 11 JJz (^ 0 5 *Hj ) 11 £00 £ A(M), 


VX|| 1)OO <||VX(0)|| liOO +A 0 f 

Jo 


IVAl I 


<||VX(0)|| ioo +A 0 


f 

Jo 


(7.95) 
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\\dt : y(a 0 ,aij)\\ LOO + Ve\\d tt yVaij\\ Lao < A(M). (7.96) 

Note that these coefficients are lack of uniform regularity with respect to normal variables. To 
get the estimates for solutions of (17.87P , the authors in [33] use the paradifferential calculous to 
prove Theorem 10.6 in [33], which yields that there exists 70 depending only on M such that 
for 7 > 70 , the solution of (17.871) satisfies the estimate 


||fi a || 2 1 < y/eMM) / \n 

Hf L 2 Jo 


a 1 2 
0 


By the Minkowski inequality and the one-dimensional Sobolev embedding H 4 c L 4 

|2 / iioai|2 / ilrAQ:il2 llrAQll2 


||o c 


IL^ 2 < 


I L^ t < ll^° 


, 1 = ||«° 

L*Hf 


HfL 2 


one then has 

\\to a \\l*L> < V^A(M) [ 

T Jo 

This and (17.861) show that 

IK||| 4l2 <a(m)V£ [ T \z a u\l 

T Jo 

Thus (17.841) follows from (17.831) and Cauchy’s inequality. 


(7.97) 


(7.98) 


□ 


8. L°° ESTIMATES 

In order to close the a priori estimates, we shall now estimate the L°° norms contained in A^, 
||9 z u|| Y m +2 and y/s \\d zz v\\ Loo , and /q \/^ ll^^lli 00 that was used in Proposition 17.51 We will 
prove that they can be bounded in terms of the quantities in the left hand side of the estimates 
of Propositions 17.31 and 16.21 which will then be shown to be in L°° in time. 

The key point is to use again the quantity S n , and one has the following: 

Lemma 8.1. For any k E N: 


\\d z v\\ Y k — A |h| Y fc+i )(l|5n|| Yfc + N| Yfc+ i). ( 8 . 1 ) 


Also, 

|| &zz v \\l 00 ^ A 0 f 

Vi II d z S n 

IIl°° ll^ n Hl,oo 1 1 ^ 1 12,00 j 

1 ( 8 - 2 ) 

and 

V^ll^v|| 1|0 o < A 0 ( 

V£\\d z s n \ 

1,00 + 2,00 + IMl3,oo) 

(8.3) 


Proof. This follows again from (17.141) and (17.161) . □ 

As a consequence of Lemma [8TT1 it suffices to estimate S n . Before proceeding further, one 
has the following observation. 

Lemma 8.2. For any smooth cut-off function x such that x = 0 in a vicinity of dQ = {z = 
0, —b}, we have for m > k + 3/2: 

llx/lliyfc’ 00 ~ \\f\\m ■ (8-4) 

Proof. It follows directly by the Sobolev embedding H k (Q) C L°°(fl) for k > 3/2 and the 
fact that away from the boundary the conormal Sobolev norm ||-|| fc is equivalent to the H k 
norm. □ 

Hence, one needs only to estimate y5 n for x( z ) compactly supported and such that x = 1 hi 
a vicinity of dkl. In this step, as in [ 32l 133] . it is convenient to use a coordinate system where 
the Laplacian has the simplest expression. We shall thus use a normal geodesic coordinate 
system in the vicinity of dfl. Note that this coordinate system has not been used before since 
it requires more regularity for the boundary: to get an H m (or C m ) coordinate system, the 
boundary needs to be Ff m+1 (or C m+l ). Nevertheless, at this step, this is not a problem since 
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one needs only to estimate a low number of derivatives of the velocity, say y + 2, while the 
boundary is H m and m can be as large as needed. We shall choose the cut-off function y in 
order to get a well defined coordinate system in the vicinity of <9H. For sake of brevity, we 
consider only the estimates near {z = 0}, and the estimates near {z = —b} may follow in the 
same way and a bit simpler. 

Define a different parametrization of the vicinity of {z = 0} by 
^(t, •) : 0 = M 2 x (—oo, 0) —> O(t) 

^• z ^{h(t,y)) +zn '’ {t ’ v '>' (8 ' 5) 

where n b = N b /|N b | is the unit exterior normal with N b = (— d\h, — d^h, 1). Note that \I/(t, ■) 
is a diffeomorphism from R 2 x (—5, 0) to a vicinity of <90(t) for some 6 which depends only on 
Co, and for every t € [0 ,T £,a ] thanks to (16.11) . By this parametrization, the induced Riemannian 
metric has the block structure 

5). (8.6) 

Hence, the Laplacian in this coordinate system reads: 

A gf = d zz f + Y z (ln \g\)d z f + Agf, (8.7) 

where \g\ denotes the determinant of the matrix g and is defined by 

Ag/ = —nr Y dyiffimfidyif)- ( 8 - 8 ) 

\9 \ 2 l<i,j<2 

To use this coordinate system, one shall first localize the equation for S^v in a vicinity of 
{z = 0}. Set 

S* = xC z)S*v, (8.9) 

where X {z) € [0,1] is smooth compactly supported near {z = 0} and takes the value 1 in a 
vicinity of {z = 0}. (J7.19I) yields that 

df S x + (v ■ \7^)S X - eA Lp S x = F S x in fi, (8.10) 

where 

Fgx = F x + F v (8.11) 

with 

F x = ( V z d zX ) S v v - eV^x ■ - eA^ x S v v, (8.12) 

F v = ~ X {D*) 2 q - | ((V^n) 2 + ((V V) 2 ) • (8.13) 

Next, define implicitly S in Q(t) by S(t,$(t,y, z)) = S x (t,y,z). Then (18.101) yields 

d t S + u-\/S - eAS = Fsxit^it,-)- 1 ) in fi(t). (8.14) 

Finally, define S ^ in H by S^(t,y,z) = z)) = S x (t, 4>(i,-) _1 o T). It then follows 

from (8.14) and (18.71) that 

d t S* + w -VS* - e(8 zz S* + ^d z (\n\g\)d z S* + A~ g S*) = Fsxit,®- 1 in Q, (8.15) 

where the vector field w is given by 

w = x(VT)” 1 (u{t, T) - d t T) = x(Vf)" 1 (v(t, 4> _1 o ) - dt*) ■ (8.16) 

Note that S ^ is compactly supported in z in a vicinity of {z = 0}. The function X (z) is a 
function with a slightly larger support such that X S^ = . The introduction of this function 
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allows to have w also supported in a vicinity of {z = 0}. Note that on {z = 0}, w 3 = 0. Indeed, 
since 

/ 1 0 n ;\ 

W(t,y,0)= 0 In* , (8.17) 

\ d\h 82 k n| J 

thus for any Y € M 3 , 

((y^(t,y,0))~ 1 Y) 3 = Y-n b . (8.18) 

Hence, (I8.16|) implies that on {z = 0}, 

w 3 = vn b - d t hn b 3 = • N fe - d t h) = 0 (8.19) 

thanks to the kinematic boundary condition. 

Now, set 

s nfay> z ) = n b (t,y)S*(t,y,z)n b (t,y) (8.20) 

with n fc = Id — n b <g>n b . Note that and n b are independent of z. Moreover, since the equation 
(18.1511 is compactly supported in z in a vicinity of {z = 0 }, this yields that S® solves 



d t Sn +w-VSn - e{d zz +-d z (ln\g\)d z )Sn = F* in {z < 0}, 

( 8 . 21 ) 

where 


Ejf = U b F S xn b + F *’ 1 + F *' 2 

( 8 . 22 ) 

with 


F *’ 1 = (( d t + w y • V y )n b )S*n b + n b S*(dt + w y ■ V y )n b , 

(8.23) 


F ^' 2 = -£U b (A~ g S qi )n b . 

(8.24) 


On {z = 0}, = S n = 0. Furthermore, it is convenient to eliminate the term ed z \n\g\d z in 

the equation (18*21 p . Set 

P=\9\*S%, (8.25) 

then 

dtp + w ■ Vp — ed zz p = TL in {z < 0 }, (8.26) 

where 

7i = \g\i (F* + F g ) with F g = p\g\~^ (d t + w • V - ed zz ) \g\F (8.27) 

Trivially, on {z = 0}, p = 0. 


8.1. Estimate of ||3 2 ,u|| Y m +2 . We now establish the first L°° estimate. Note that it is equivalent 
to estimate S n or or p. Indeed, by the definition (18.2011 . using the chain rule and the fact 
that Z is tangent to {z = 0}, one has 


I'-’n 11Y fc — ^ 


CO 


l/llyfe + I 


n b S*v n b 


Y k 


Since |n — n 6 | + |n — n 6 | = O(z) in the vicinity of {z = 0}, thus 


H'S’n || y* “ ^ (co ’ ^ ¥fe+1 ^ (11 11 + IIC || Y fc +! ) ■ (8.28) 

Similar arguments show that 

11‘S'n11 Y fc — A l^'lY fc + 1 ^ (ll^n IIy* + IMIy*^ 1 ) • (8.29) 


On the other hand, it is easy to see that it is equivalent to estimate p or . By (|8.1H . it thus 
suffices to estimate p. 

Set 

—) H x ^+7 + ll^^ll y ^+2 + IMIxif +7 + ll^ 7t ’llx^ :+6 


Aoo(i) = A 


(8.30) 
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Proposition 8.3. For m > 14, it holds that 


\\p( t )\\ Y I T+ 2 ~ IIp( 0 )II y ?+ 2 + 


f A 00(1 

J 0 


+ £ || VS* n ) 


Proof. Apply Z a . for a € N 1+3 with |a| = k < Pf + 2, to (18.2611 to obtain 
d t Z a p + w • VZ a p - ed zz Z a p = Z a PL + C s in {z < 0}, 

where 


Cs=Cl + C 2 




with 


Cs = [Z a , Wy} ■ V y p + [Z a , u; 3 ] • d z p := C Sy + C S z and C 2 S = -e [Z a , d zz \ p. 
The maximum principle on (I8.32|) yields that, since Z a p = 0 on {z = 0}, 


(8.31) 

(8.32) 

(8.33) 

(8.34) 


\Z a p\\ L oo < \\Z a p(0)\\L~ + f (\\Z a n L °° + \\Cs\\ L -) • (8-35) 

Jo 


The right hand side of (18.351) can be estimated as follows. For the commutator Cg, the direct 
estimates yield 

II^Sj/ll^oo ll u; j/llY fc IIpIIy* ■ (8.36) 

To estimate Cs z , by expanding the commutator and using (13.111) . one needs to estimate terms 
of the form 

| zP W3 d z z^p 

with (3 + 7 < a and |y| < |a| — 1. Since w 3 = 0 on {z = 0}, so 

Z^W 3 


Z ls w 3 d z Z' 1 p 


z{z + b ) 


Z 3 Z~<p 


< 

rs_/ 


||9 2 U;3|| Y fc ||p|| Y fc • 


For the commutator C|, using (13.111) repeatedly leads to 

- [z a ,d zz ] P = d z ([z“, d z ] P ) + [z a ,d z \ d zP 


(8.37) 

(8.38) 


^2 9 Z 

\/3\<k-l 


+ ^2 c pd z (Z l3 d z p) 

m<k-i 


Hence, 


= ^2 c /3 >Z l3 ’d z p+ ^2 cp'Z P 'd zz p. 

| 0 '|<fc-l \(3'\<k-l 

l^illioo ~ e \\9zP\\ Y k-i + e \\9 zz p\\ Y k-i ■ 


Note that using the equation (18.261) implies 

£ \\d zz p\\ Y k-i < ||9j/? + W-Vp 7^llyfc-i ^ ||p||Y fc d - 11^ ■ Vp|| Y fc_i l|7^ || Y fc-i 

Recall from (17.11) that 


As (18.361) and (18.371) . 


-1 < eA ( —, \h\ Yk +i ) ||V ^|| Yfc _i < Aoo. 


\W ■ Vp|| Y fc-l < (Ikylly* + ||3 Z W 3 || Yfc ) ||p|| Y fc 


Thus, 


C 


S11 ^00 


(l + Hwj/IIy*; + 11 <9z ^ 3 11 Y fc ) IIpIIy* "b II^IIy*- 1 + Aqo. 


(8.39) 

(8.40) 

(8.41) 

(8.42) 

(8.43) 


Consequently, in light of the estimates (18.361) . (18.371) and (18.431) into ([8.351) . one obtains 


WpWy* ~ IIp(O) II Y fc + / II^IIy* + (■*■ + II w 2 /ll Y fc + ll 5 *“3 
Jo 

Now we estimate PL. Note first that 


+ A c 


\g\ 4 Fg Yfe < A , ||p|| Y fc + IMIy* + Wy*+3^ • 


(8.44) 
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Next, it follows from (18.231) and (18.241) that 

Idlin ' 11 < A I^Iy^+2 + IMIy* + II Vu IIy^ 


and 


4 17^,2 


\g\*F* 


Y k 


< eA I —, |/i| ¥ fc+3 ) ||Vu|| Y * + 2 • 


Co 


Using (18.121) and (18.131) . the fact that F x is supported away from {z = 0} and Lemma [8721 one 
gets 


| ff |4lI 6 FV 


Y k 


< A 


\h\ Y k +2 + ||p|lx fc + 4 


and 


|<7|3n^n b || Yfc < A Q-, |h| ¥ , +2 + ||Vu|| ¥ fc^ (1 + \\V 2 q\\ Yk ) . 
Recalling (8.27), (8.22) and (8.11), and collecting these estimates, one arrives at 
II'HHy*, < A , |/i| Y fc +3 + |M| X fc + 4 + ||Vu|| Y fc^ (l + £ ||Vu|| ¥ fc +2 + || Vug'll Yfc ) . 
Recall from (14.151) that 


(8.45) 


II^^IIyTT +2 — ^ | ^-| y Z? 2+6 +4 + |h| x m+6 +4 + || Vu || ^ m+6 +3 + 11 11 ^ rri+S +4 + 11 Vu 11 ^ rn +6 +3 ^ 

< A , |h| Y m +5 + |h| x ™, +7 + ||Vu|| y ™, +2 + ||v|| x ^+7 + l|Vu|| x ^:+6^ (8.46) 

if m > 10. Hence, (18.441) implies 

HpII Y ^+2 < ||/?o|| ¥ ^+2 + j A , |^l x ^+7 + ||Vu|| Y m +2 + ||u]| x m +7 + 11 Vu 11 x ™ +6 (8.47) 

x (l + e||Vu|| Y - +4 ) . 

Then the desired estimates (8.29) follows for m > 14 so that y + 4 + l<m — 2. □ 


8.2. Estimate of \f£ \\d zz v\\ Loo . The next L°° estimate is the only place where one needs to 
use the compatibility condition on the initial data. As in the previous subsection, by (18.21) . one 
can reduce the problem to the estimate of \fe ||5 z p|| ioo . 


Proposition 8.4. Assume that the initial data satisfies the compatibility condition S’ n (0) = 0 
on {z = 0, —b}. Then it holds that for m > 6, 


Vs\\d z p{t)\\ Loa < Ve\\d z po\\ L , 


f l Aqc 

Jo \Jt-T 


(8.48) 


Proof. The proof follows the spirit of the proof of Proposition 9.8 in [33|. Recall that p satisfies 
(18.261) in {z < 0} with p = 0 on {z = 0}. Note that one can not apply \fed z to (18.261) and then 
use the maximum principle due to boundary condition. We shall use a precise description of 
the solution of (18.261) . Indeed, one can use the one-dimensional heat kernel of {z < 0}: 


G(t, z, z') = 


\J An et 


e 4et — e 4et 


(8.49) 


to write that 


r+o o 

y/ed z p(t,y,z) = / \fed z G(t, z, z')po(y, z')dz' 

Jo 

+ f \fed z G{t — t,z,z') ( TL(T,y, z') — w ■ Vp) dz'dr,. 

Jo 


(8.50) 
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Since po = 0 on {z = 0}, thanks to the compatibility condition, one can integrate by parts the 


first term to obtain 

V£\\dzP(t)\\ L oo<Ve\\d*Po\\ L °°+f /t~— (I|7*IIl°° + Ik ■ Vp|| L oo) ■ (8.51) 

Jo ~ T 

Next, it follows from (I8.45|> with k = 0 that 

\\H\\ LOO < Aoo (1 + ||V 2 g|| Loo ) < Aoc. (8.52) 

On the other hand, as (18.421) . 

Ik • Vpllioc < (Ikyllioc + |k^3||l,oc) ||p||l,oo ^ ( 8 - 53 ) 

Consequently, plugging (18.521) and (18.531) into (18.511) yields that 

Ve\\dzP(t)\\ Loo < Ve\\d g po\\ Lao + f (8-54) 

Jo y/t-T 

which completes the proof. □ 


8.3. Estimate of yfz ||<9 zz u|| i oo- The last L°° estimate is the one that was used in Proposition 
17.51 By (18.3D . one can again reduce the problem to the estimate of y/e ||<9 z p|li 00 . 

Proposition 8.5. For m > 6, it holds that: 

[ >/e\\dzP\\i t00 £y/i\\po\\ loo + t f -}== ■ (8-55) 

Jo Jo V t — T 

Proof. We will use a different argument from the proof of Lemma 9.9 in |33| . A direct use of 
the Duhamel formula (18.SOD yields 

Ve|IMi,oo^^llA)lli,oo + ^ (H^Hi.oo + Ik ~ Vp|| 1)00 ) ■ (8.56) 

Next, it follows as the previous arguments that 

Ikllgoo < Aoo (l + ||V 2 g|| l oo ) < Aoo (8.57) 

and 

Ik • Vp|| 1)QO < Aoo (lkyll 2 ,oo + lk^3|| 2 ,oo) IIpII 2 ,oo < Aoo- (8.58) 

Consequently, 

^IIAplli.oo S ^ IIpoIU.oo + J‘ ^=- (8-59) 

Integration in time yields 

[ Vz\\dzp\\ lt00 <Vi\\po\\ loo + t f -}== . (8.60) 

Jo Jo y/t-T 

which completes the proof of the proposition. □ 

9. Proof of Theorem nm 

Collecting the estimates obtained in Sections [6][8l one can prove Theorem 12.11 in the similar 
way as that for Theorem 1.1 of [33] with slight modifications. 

Recall J\f(T) and Q(T). For two parameters R and co to be chosen 1/cq R , define 

Tf’ S = sup It € [0,1] | M{t) < R, \h{t)\ 2 < —, d z <p(t) > ^ (9.1) 

T>0 l Co 2 

and g — dfq(t ) > on {z = 0}, Vi € [0,T].| . 

Proposition 18.51 yields 



(9.2) 
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This allows one to use Proposition 17.51 which together with Propositions 16.2116.3117.3118.31 and 
18.41 implies that, by a suitable linear combination, 


Af(t) < A f ^-,Af(0)^ + MR) ( T * + VW(T) + W(T)s) , 


where 


co 

T 


W(T) := [ \d?h\ 2 0 + a \d?h\l + \\d?v\\ 2 + \\d z v\ 
Jo 

It follows from the Cauchy-Schwarz inequality that 

W(T) < A7(T)T^ < RTF 


Hence, one deduces from (19.31) that 


A f(t) < A (^—,AA(0)J + A(i?)T3. 

On the other hand, since J\f(T) involves time derivatives, one gets easily that 


1 


IMi)l2,oo<IM0)| 2lO o + A(i?)T, 

d z <p(t) > dMO) - A (R)T 

and 

9 — dfq(t) > g — dfq(0) — A(R)T. 

Consequently, one can choose Co so that |/i(0)| 2oo < 277 an d then R = 2A 
there exists T* which depends only on R so that for T < Min (T*,T* ,<5 ), 



(9.3) 


(9.4) 

(9.5) 

(9.6) 

(9.7) 

(9.8) 
then 


A f{t) < d x <p(t) > \h(t)\ 2 oo < and g - dfq(t) > ^ on {z = 0}, Vt € [0,T], 

This yields Tf' S > T* by the definition (19.11) and also the estimate (|2.6D . The proof of Theorem 
o is thus completed. □ 


Appendix A. Sobolev conormal spaces 

We recall the Sobolev conormal spaces X m and Y m from m- 

Lemma A.l. The following product and commutator estimates hold, 
(i) For |or| + |/3| = k > 0; 


k > 1: 

z a /zfg < ||/|| x » |I»II y 5 + II/II y j llsllx* ■ 

(A.l) 

\\[Z a J]g\\ < WZfW^ ||^|| yV + \\Zf\\ Yh -i || 5 ||x*-i 

(A.2) 

k> 2, 

define the symmetric commutator 



[Z a J,g} = Z a (fg)-Z a fg-fZ a g. 

(A.3) 

IP“ 

f,g] II < \\zf\\*~* \\Zg\\ Y ^ + \\zf\\ Y k^ \\Zg\\ xk - 2 . 

(A.4) 


Then 


Proof. The product estimate (Bud follows by controlling the product with the lower order 
derivative term in L°° and the higher order derivative term in L 2 . To prove the commutator 
estimate ED , one uses the Leibnitz formula to expand 


[Z a ,f]g= ]T C^fZHg. 

/3+7=q: 

/ 3#0 

Since /3 / 0, one can write Z$ = Z^~^'Z&' with \(5'\ = 1. Then (1A.1D yields 
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This proves (1A.2I) . The commutator estimate (1A.4D can be proved in the same way. 
We shall also use the Sobolev tangential spaces defined by 

HL n m = {/ e L 2 (n), \\f\\ HLn = ||A S /|| L2 < 00 } , s € R, 

s 

where A s is the tangential Fourier multiplier by (l + |£| 2 ) 2 , £ G R 2 . Note that 


□ 


Hi 


< 


k for s < k, k £ N. 


Lemma A.2. The following anisotropic Sobolev embedding and trace estimates hold, 
(i) For si + s 2 > 2, s 3 + s 4 > 2: 


LOc < \\d z f\\l si ll/ll 


1 

2 

H s 3 
11 tan 


(ii) For si + S 2 = S 3 + S 4 = 2s: 


I/I 


s < \\d z f\\ 2 H n 


1 

2 

tt s 2 
11 tan 


1 

2 

11 tan 


1 

2 

tt s A 
11 tan 


(A.5) 

(A. 6 ) 


Proof. We need to modify the proof of Proposition 2.2 in [33] since our domain here is of finite 
depth. To get the anisotropic Sobolev embedding estimate (1A.5[) . one first notes that 

/(£,z) < f(£,z') + (/ , 2 d ~f(C,x 3 ) f(£,x 3 ) dx 3 ^J . 

Integrating the inequality above with respect to z' € (—b, 0) yields 


f(€,z) ^ j ^ f(£,z') d z ' + yj b dzf(f,,x 3 ) f(£,x 3 ) 


dx 3 


Hence, it follows from the Cauchy-Schwarz inequality and the fact that si + S 2 > 2, S 3 + S 4 > 2 
that 


L<x < sup 


zG(—6,0) ./R 2 
r -0 


/(£>*) 


df 


< 

r^j 


[ [ \^ z ' 

J R 2 J-b 


dfdz' + b 


i 2 \J-b 


dzf{f,,x 3 ) /(£,® 3 ) 


dx 3 df 


I —b ■ 


(i + ieir 3+S4 f&z 


dzdf 


dfdz j +(f (1 + |£| ) S1+S2 ( dj(f,z) f(£,z) 

J ]R 2 */ b j 

< ||A S3 /||i ||A S4 /||i + ||^A Sl /||i ||A S2 /||i . 

To prove the trace estimate (1A.6[) . since si + S 2 = 2s, one may write 

\f(-iti)\ 2 Hs = \f{-iz')\ 2 Hs +[ [ 2d z A s f(z,y)A s f(z,y)dzdy (A.7) 

J R 2 Jz' 

= \f{-,z')\ 2 H3 + [ f 2d z A sl f(z,y)A S2 f(z,y)dzdy 
J R 2 Jz' 

Integrating the equality above with respect to z' G (—b, 0) and using the Cauchy-Schwarz 
inequality give the desired estimate. □ 

Following similar arguments, we also have the following Poincare inequality. 

Lemma A. 3. It holds that 

ll/ll <l/l„ + l|0,/||. (A.8) 

Proof. The proof of the estimate (lA.fij) with si = S 2 = s = 0 also leads to 

ll/ll <l/lo + ll^/ll § ll/ll*- 

Then the Poincare inequality (IA. 8 D follows by Cauchy’s inequality. □ 
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We also recall the classical product and commutator estimates in R 2 : 

Lemma A. 4. The followings hold. 

|A s (/ff )| L 2 < |/| £ oo |A s < 7 | L 2 + \g\ Loa |A s /| l2 for s > 0, (A.9) 

I [A 5 ,/] V y g\ L2 < \V y f\ LOB \A s g\ L2 + \V y g\ Loo \A s f\ L2 for s > 0 , (A. 10 ) 

l/s|l 0 l/li,00 Iff I i and l/ 5 l_i Z I/ll,00 bl_i • (A.ll) 

Proof. These estimates (1 A.9jl and (1A.10D are classical, see |2S] for example. Note that 

I/5I0 ~ bio and I/5I1 ;$ l/li,00 bli, 

the estimate (jA.Ill) follows by the theory of interpolation and duality. □ 

Note that Lemma lA.il also holds on R 2 , while we also need the following for half regularities. 
Lemma A. 5. For |a| + |/3| = k > 0: 

1 1 


Z a fZPg 


< 


I/Ixm bLfe+i + 1 / 1 * 


s ~ '^'y3 +1 j y^ +1 


s = — - 


2 ’ 2 


(A.12) 


Proof. The estimate (IA.12D follows by using (1A.11I) to control the product with the higher order 
derivative term in H s and the lower order derivative term in W 1,oa . □ 

Appendix B. Poisson extension 

We recall the extension r] of h onto {z < 0} defined by (11.61) with parameter A > 0 in the 
following form 

(B.l) 


V(y,z) = (1 + -J C(y,z) with = exp(A|£|z)/i(£). 

We hrst verify that <p dehned by (11.51) is a diffeomorphism. 

Proposition B.l. Assume that h$ £ H s ({z = 0}),s > 5/2 and ho > — b. Then there exists 
sufficiently small A > 0 such that 


1 


1 


d z tp 0 > - 1 + -ho > 0 in Q 


(B.2) 


[ d 2 z ((£,x 3 ) 

Jz’ 


dx 3 . 


(B.3) 


Proof. Note that 

dz({£,z) < d z ((£,z') 

Integrating the inequality above with respect to z' £ (—6,0), one can deduce 

6||<9 2 CIIloc < 6 sup [ d z ( (£,z) df < [ [ d z ({f.,z) dz + [ [ d%C(€,z) 

ze[-b, o]^r2 j_ b J-bJmf 

For s > 1, it then follows from the Cauchy-Schwarz inequality and the definition (IB. II) that 


dz. 


f [ hMz) 

■J~b J R? 


dz < 


r0 [ (i + bi) 2s dMzf 


dz 


= A 


= A 


(i + birw ho 


(1 + blHei 2 HO 


J exp [2A\£\z)dz 
2 1 — exp (2A\£\z) 


2 Ab| 


dz 


and 


f [ d z aoz) 

J-b J K? 


dz< 


(i + bir d z aoz) 


' —b ■ 


dz 


A? \h\ 


S +l 
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= A 2 


(1 + ICI) 2s |C| 4 HO 2 f° exp (2A\l\z)dz ) 2 < \h\ s+ s 

! J-b J 


We thus deduce 

\\d z (\\ Laa <A^\h\ s+l , 

Now we prove (|B.2[) . It follows from the definitions (11.51) and (IB. II) that 


(B.4) 


d z <po — 1 + ^Co + ^1 + d z (o. (B.5) 

By (IB. 41) . this yields that for s > 

dzP o > 1 + -^ho + — (Co — ho) + ^1 + — ^ 9 2 Co 

> 1 + — ^2 + — ^ H^zCo||x,oo > 1 + —ho — A 2 ^2 + — ^ |ho|| > — ^1 + —ho^ 

if A has been chosen sufficiently small. □ 

We also have the following well-known estimates for 17 . 

Lemma B.2. For s € M: 

Ihllj/a < |h| s _i . (B. 6 ) 

For k € N: 

IMI w k ’°° ~ l^lfc.oo 1 (^-7) 

Proof. One deduces in the same way as Proposition 3.1 in [33 that 

IIveil H* £ \h\s+± » s € M and h\\ W k,oo < \h\ koo , k € N. 

Then the estimates (IB. 6 I) (IB.71) follow by noting that rj = (1 + |)C for z € [—6, 0]. □ 


Appendix C. Some geometric estimates 

We recall that the control of quantities like f\ 2 dVt yields a control of the standard H 1 

norm of /. 

Lemma C.l. Assume that d z (p > Co and ||Vc/ 2 || L oo < ^ for some Co > 0, then 

||V/|| 2 < A 0 / |V^/| 2 dV t . (C.l) 

Jn 

Proof. We refer to Lemma 2.8 in [33] • Cl 

We also need the Korn type inequality to control the energy dissipation term. 

Lemma C.2. Assume that d z (p > Co and ||V</?|| ioo + ||V 2 <y?|| Loo < ^ for some Co > 0, then 

||Vu|| 2 <A 0 (^j \S^v\ 2 dV t + \\v\\A . (C.2) 

Proof. We refer to Proposition 2.9 in [33]. □ 

Finally, we will also need the following H _1 / 2 boundary estimates for functions satisfying 
v € L 2 and V v • v € L 2 . 

Lemma C.3. If ||V</?|| LOO < y- for some cq > 0, then 
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Proof. We adapt the proof of Lemma 3.3 in m- We will only prove the result on {z = 0}. Let 
if € H 1 / 2 , and let ip € iL 1 (n) be a bounded extension. Then 


[ ipv-N = f 

J z =0 J !T2 


ipv • N = / V* ■ ('ipv)dVt = / (■ V + ipvv -v) dVt 


< A 0 


ip 


||V*"u|| + 


Vt/> 


^||) < A 0 \ip\i (||u|| + ||• v\ 


Then the estimate ()C.3H follows from this inequality above by taking the supremum over all if 
so that |t/’|i < 1. □ 
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